APERIODIC OSCILLATORY ASYMPTOTIC BEHAVIOR 
FOR SOME BIANCHI SPACETIMES 



FRANgOIS BEGUIN 

Abstract. We study the asymptotic behavior of vacuum Bianchi type A spacetimes close to their singu- 
larity. It has been conjectured that this behavior is driven by a certain circle map, called the Kasner map. 
As a step towards this conjecture, we prove that some orbits of the Kasner map do indeed attract some 
solutions of the system of ODEs which describes the behavior of vacuum Bianchi type A spacetimes. The 
orbits of the Kasner map for which we can prove such a result are those which are not periodic and do 
not accumulate on any periodic orbit. This shows the existence of Bianchi spacetimes with an aperiodic 
oscillatory asymptotic behavior. 



1. Introduction 

1.1. Bianchi spacetimes, tlie Wainwriglit-Hsu vector field and tlie Kasner map. A Bianchi space- 
time is a cosmological spacetime which is spatially homogeneous. More precisely, it is a spacetime 
(M, g) such that M is diffeomorphic to the product G x / of a simply connected three-dimensional Lie 
group G and an interval / of the real line, and g = ht — dt^ where ht is a left-invariant riemannian 
metric on G x {t} ~ G for every t & I. A vacuum Bianchi spacetime is a Bianchi spacetime (M, g) 
satisfying the vacuum Einstein equation Ricg = 0. A type A Bianchi spacetime is a Bianchi spacetime 
for which the corresponding three dimensional Lie group G is unimodular. Note that the Lie group G 
can be assumed to be simply connected without loss of generality. 

A Bianchi spacetime can be described as a one-parameter family {ht)tei left-invariant riemannian 
metrics on a three-dimensional Lie group G. The space of left-invariant riemannian metrics on a given 
Lie group is finite-dimensional. Therefore, when restricted to the context of Bianchi spacetimes, the 
vacuum Einstein equation becomes an ODE on a finite dimensional phase space. A classical and conve- 
nient way to write this ODE is to use the so-called Wainwright-Hsu variables (Si, S2, S3, Ni,N2, N3). 
Very roughly, Ni,N2, are the structure constants of the Lie algebra of the Lie group G x {t} in a 
certain basis, and Si, S2, S3 are the components of the normahzed traceless second fundamental form 
of G X {t} in the same basis. The corresponding phase space is the four-dimensional manifold 

B := {(Si,S2,S3,iVi,Ar2,iV3) eM^ | Si + S2 + S3 = , = 0} 

where 

= 6 - (S? + si + si) + ^(iVf + iV| + iV|) - (N1N2 + iViAf3 + N2N3). 
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The vacuum Einstein equation is equivalent to the system of ODEs 



(1) 



where 



and 





= (2-Q)Si-i?i 




= (2-g)S2-i?2 




= (2-(?)S3-i?3 


N[ 


= -(g + 2Si)iVi 




= -(g + 2S2)Ar2 




= -(g + 2S3)iV3 


Q = 





' Ri = I (27V2 - 7V| - iV| + 2iV2Ar3 - N1N3 - N1N2) 

R2 = l{2Ni-Ni-N'f + 2N3Ni-N2Ni-N2N3) 

_ R3 = l{2Nl-N^-Nl + 2NiN2-N3N2-N3Ni) 
(see e.g. [3] for the construction of the Wainwright-Hsu variables and the expression of the Einstein 
equation in these variables). In other words, vacuum Bianchi type A spacetimes can be seen as the 
solutions of the system of ODEs (1) on the four-dimensional manifold B. We will call Wainwright-Hsu 
vector field, and denote by Xq, the vector field on B associated to the system of ODEs (1). We will 
denote by the time t map of the flow of Xg. To measure the distances on the phase space B, we will 
use the riemannian meti-ic /i = (dEi)^ + {dT,2)^ + (d'^s)^ + (dNi)^ + {dN2)'^ + {dN^f. 

Remark 1.1. We have chosen the "anti-physical" time orientation. With this convention, Bianchi space- 
times are future incomplete, but not necessarily past incomplete. The main reason for this choice is that 
we want the so-called mixmaster attractor (see below) to be an attractor, rather than a repellor. 



The phase space B admits a natural stratification, which is invariant under the flow of the Wainwright- 
Hsu vector field Xq. There are six strata denoted by Bl,BJl,B^fl^,B^fJl^,B^fJ^,Bl^. These strata 

correspond to the different possible signs for the variables A'^i, N2, A^a- They also correspond to the 
different (isomorphism class of) simply connected unimodular three-dimensional Lie groups. The orbits 
of Xb contained in Bi are called type I orbits, the orbits contained in Hjj are called type II orbits, etc. 
The behavior of the type I, II, VIq and VIIq orbits of Xq is very well understood. On the contrary, the 
behavior of the type Vlll and IX orbits (which are the generic orbits in B) is, at best, conjectural. In 
order to simplify the discussion we will focus our attention on the subset B^ of B where the coordinates 
Ni,N2, N3 are non-negative: 

B+ := {(Si, S2, E3, m,N2, iV3) G B I TVi > , 7V2 > , iV3 > 0}. 

Observe that this subset is invariant under the flow of X^. 

The stratum Bj is more frequently denoted by /C. This is a euclidean circle in ;B C M^, called the 
Kasner circle. It is made of the points of B where Ni = N2 = Ns = 0. This corresponds to the case 
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where the Lie group G is abelian {i.e. G = M?). 

(2) /C = {(Si,S2,S3,Afi,iV2,A^3) G^l A^i =A^2 = iV3 = 0} 

= {(Si, S2, S3, iVi, iV2, A^a) G I A^i = iV2 = iVa = , 

Si + S2 + S3 = , S? + S2 + S2 = 6} 

The vector field vanishes on JC. Therefore every type I orbit is a fixed point. There are three 
points on /C, called the special points or the Taub points, which play a very important role in every 
attempt to understand the behavior of the solution of the Wainwright-Hsu equations; these are the points 
for which (Si, S2, S3) is equal respectively to (2, —1, —1), (—1, 2, —1) and (—1, —1, 2); these points 
are usually denoted by Ti,T2,T^. For every point p ^ JC, the derivative DXb{p) has four distinct 
eigendirections: the direction TpK. of the Kanser circle corresponds to a zero eigenvalue (since Xq 
vanishes on /C), and the directions M.g^(p), M.Tj^(p), R.g^(p) are respectively associated to the 
eigenvalues -(2 + 2Si), -(2 + 2S2), -(2 + 2S3). When p is one of the three special points Ti, r2, Ta, 
two of the three eigenvalues —(2 + 2Si), —(2 + 2S2), —(2 + 2S3) vanish; hence, the derivative of 
DXsip) has a triple zero eigenvalue. If p is not one of the three special points Ti,r2,r3, then the 
eigenvalues —(2 + 2Si), — (2 + 2S2), —(2 + 2S3) are pairwise distinct, two of them are negative, and 
the other one is positive; in other words, the derivative DXq[p) has two distinct negative eigenvalue, a 
multiplicity one eigenvalue, and a positive eigenvalue. 

The stratum Bjj is two-dimensional. It is made of the points of B for which exactly two of the A?j's 
vanish. This corresponds to the case where the Lie group G is isomorphic to the Heisenberg group. 

BiU% = {(Si,S2,Sa,iVi,A^2,A^3)eH| A^i = A^2 = 0oriVi = Ara = 0or7V2 = 7Va = 0} 

One can easily check that Bi U Bn is the union of three ellipsoids which intersect along the Kasner 
circle. The vector field Xg commutes with the transformation (A'^i, A'^2, N^) (— A'^i, — A^2, —-^3), so 
we can work in Bn H B^ to simplify the discussion. The intersection Bn n B^ is the union of three 
disjoint hemi-ellipsoids bounded by the Kasner circle. Type II orbits can be calculated explicitly (see for 
example [3]). It appears that every type II is a heteroclinic orbit connecting a point of /C \ {Ti, T2, Ta} 
to another point of /C \ {Ti, T2, Ta}. Moreover, for every point p € /C \ {Ti, r2, T3}, there is exactly 
one type 11 orbit in B^ which "takes off dXp (this orbit is tangent to the eigendirection associated to the 
positive eigenvalue of DXb{p)) and two type II orbits in B^ which "land" at p (these orbits are tangent 
to the eigendirections associated to the negative eigenvalues of DXb{p)). 

The Kasner map / : /C ^ /C is defined as follows. Consider a point p G /C. If p is one of the 
three special points Ti, T2, T3, then f{p) = p. Otherwise, one considers the unique type II orbit which 
"springs up" at p; this orbit connects p to another point G /C; this point p' is by definition the image of 
p under /. For every p G /C \ {Ti , T2, T3}, we will denote by Opj(^p^ the unique type II orbit connecting 
the point p to the point f{p) in B'^. Since there are exactly two type II orbits in B^ which land at a given 
point p ^ K,, the Kasner map is two-to-one. The exact computation of the type II orbits of Xjs yields 
a nice geometric description of the Kasner map /. Consider the equilateral triangle which is tangent 
to K at the three special points Ti, T2, r3. Denote by Mi, M2, M3 the vertices of this triangle. For 
p G /C \ {Ti, r2, T3}, let Mi be the vertex of /C which is the closest to p. The line {p Mj) and intersects 
the circle K, at two points: the point p and the point f{p). See figure 1. Using this geometric description, 
it is easy to see / is C^. Moreover, the Kasner map / is non-uniformly expanding. Indeed, the norm of 
the derivative of / (calculated with respect to the metric on /C induced by the riemannian metric h) at a 
point p G /C is equal to 1 if p is one of the three special points Ti , r2 , Ta, and is strictly bigger than 1 if p 
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Figure 1 . Construction of an orbit of the Kasner map 

is not one of these three special points. It follows that, for every compact set C in /C \ {Ti, T2, T3} there 
exists a constant f c > 1 such that 1 1 1 D/ (p) 1 1 1 > z^c for every p in C. 

The strata Byi^ and Byn^ are three-dimensional. The stratum Byi^ is made of the points of B for 
which exactly one of the A'j's is equal to zero, the two others being of different signs. This corresponds 
to the case where the Lie group G is isomorphic to Isom(R^'^) = 0(1, 1) k M^. This stratum is disjoint 
from B^. The stratum ^vilo made of the points of B for which exactly one of the A^j's is equal to zero, 
the two others being of the same sign. This corresponds to the case where the Lie group G is isomorphic 
to Isom(M2) = 0(2) >< r2_ 

Finally the strata .Svni ^nd Bix are four-dimensional (these are open subsets of B). The stratum ^vill 
is made of the points of B for which the Ni 's are non-zero and do not have the same sign. It corresponds 
to the case where the Lie group G is isomorphic to the universal cover of SL(2,M). It is disjoint from 
B^. The stratum Bix is made of the points of B for which the Ni's are either all positive, or all negative. 
This corresponds to the case where the Lie group G is isomorphic to S0(3, M). 

1.2. Statement of the main results of the paper. Misner has conjectured that the dynamics of type 
IX orbits of the Wainwright-Hsu vector field Xq "is driven" by the Kasner map ([5]). The idea is that 
every type IX orbit should eventually approach the so-called mixmaster attractor A:= KVJ Bn, and then 
"follow" a heteroclinic chain q — )- Oqjt^q) — )• f{q) — )• C'j(^) j2(g) —)-.... Conversely, every heteroclinic 
chain as above should attract some type IX orbits of Xb- Moreover, "generic" type IX orbits should be 
attracted by "generic" heteroclinic chains. This would imply that the behavior of generic type IX orbits of 
Xq is determined by the behavior of generic orbits of the Kasner map. See [3] for a detailed discussion 
of various possible precise statements for this conjecture. 

In 2001, Ringstrom proved that A := K,\J B\i is indeed an global attractor: the distance from almost 
every type IX orbit of the Wainwright-Hsu vector field Xjs to A tends to as the time goes to 00 ([9], see 
also [4]). Ringstrom's theorem has important consequences, such as the divergence of the curvature in 
Bianchi spacetimes as one approaches their past (or future) singularity. Nevertheless, this theorem does 
not tell anything on the relation between the dynamics of Bianchi orbits and the dynamics the Kasner 
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map. The purpose of the present paper is to prove that "there are many points q E JC such that the 
heteroclinic chain q O^jf^q^ f{q) ^f{q),p{q) ■ ■ ■ does attract some type IX orbits of Xb". 

Definition 1.2. For g G /C, we call stable manifold ofq, and we denote by W^{q), the set of all points 
r G ;B for which there exists an increasing sequence of real numbers (tn)n>o such that 

dist(X^"(r),r(g))^-H^O 

and such that the Hausdorff distance between the piece of orbit {^^(r) ; tn <t < tn+i} and the typell 
heteroclinic orbit Ojn(^) jn+i(q) tends to as n goes to +oo. 

A subset C of /C is said to be forward-invariant if /(C) C C. It is said to be aperiodic if it does not 
contain any periodic orbit for /. Our main result can be stated as follows: 

Theorem 1.3. If q G K. is contained in a closed forward-invariant aperiodic subset of /C, then the 

intersection of the stable manifold W'^{q) with Bjx is non-empty. More precisely, W^{q) H Bjx contains 
a embedded three-dimensional closed disc. This disc can be chosen to depend continuously on q (for 
the topology on the space ofC^ embeddings of the closed unit disc ofW^ in B) when q ranges in a 
closed forward-invariant aperiodic subset ofK. 

Remark 1.4. Actually, one can prove that W^{q) fl Hjx is an injectively immersed open disc, and that 
this depends continuously on q (for the compact-open topology on the space of immersions of the 
open unit disc of in Bix) when q ranges in a closed forward-invariant aperiodic subset of /C. Note 
that (q) n Bix is never a properly embedded open disc. 

Remark 1 .5. We decided to focus on type IX orbits for sake of simplicity. Nevertheless, the analog of 
theorem 1.3 for Bym instead of Bix is also true. The proof is exactly the same: one just needs to replace 
the set B+ by the set :=Bn {Ni < 0} n {N2 > 0} n {N3 > 0}. 

Observe that the hypothesis of theorem 1.3 is satisfied by a dense set of points in /C : 

Proposition 1.6. The union of all the closed forward-invariant aperiodic subsets ofK is in dense in /C. 

Remark 1 .7. Let E be the union of all the closed forward-invariant aperiodic subsets of fC. According to 
the above proposition, E is dense in /C. Observe nevertheless that is a "small" subset of /C both from 
the topological viewpoint (it is a meager set) and from the measurable viewpoint (it has zero Lebesgue 
measure). Also observe that theorem 1.3 does not tell that W'^{q) depends continuously on q when q 
ranges in E. So, we do not know whether the union of all the stable manifolds W^{q), where q ranges 
E, is dense in B (or in an open subset of B) or not. 

The origin of the hypothesis of theorem 1.3 is purely technical: if p belongs to a closed forward- 
invariant aperiodic subset of /C, we can find a coordinate system in which the vector field Xq depends 
linearly of all the coordinates but one; this makes the estimates on the flow of Xjs near such a point much 
easier. We do not know if such Unearizing coordinate systems exists near a point of /C which is periodic 
or preperiodic under the Kasner map. Nevertheless, it should be noticed that the estimates on the flow 
of Xb that we need to construct stable manifolds are much weaker than the existence of a linearizing 
coordinate system. So, we do think that theorem 1.3 can be extended to any point g € /C such that 
the orbit of q under the Kasner map does not accumulate one any of the three special points Ti, Ta. 
Understanding the behavior of the orbits which pass arbitrary close to the three special points seems to 
be a much harder problem. 
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Remark 1.8. A consequence of theorem 1.3 (and proposition 1.6) is that the Wainwright-Hsu vec- 
tor field is, at least partially, sensitive to initial conditions: there is a dense subset E of the Kas- 

ner circle /C such that, for every point q ^ E, arbitrarily closed to q, one can find two points 
ri , r2 G Bix such that the orbits of ri and r2 will not have the same "asymptotic behavior" (for ex- 
ample, limsupdist/i (Xg(ri) , X^(r2)) > 1/10. 

t-^+oo 

As we were finishing to write the present paper, M. Georgi, J. Harterich, S. Liebscher and K. Webster 
put on arXiv a preprint in which they prove that the period three orbits of / admit non-trivial stable 
manifolds ([2]). At the end of this preprint, they claim that their techniques can be used to extend their 
result to any periodic orbit of /, and even to any orbit of / which does not accumulate one any of the 
three special points Ti, T3. Such a extension would imply our result. 

1.3. Idea of the proof and organization of the paper. Let us try to sketch the key idea of the proof 
of our main theorem 1.3. Let C be a closed forward- invariant aperiodic subset of the Kasner circle, 
and denote by C the union of C and all the type 11 orbits connecting two points of C in B^. We will 
construct a kind of Poincare section adapted to C: a hypersurface with boundary Mc which intersects 
transversally every type II orbit connecting two points of C in B~^. Theorem 1.3 will follow from the 
existence of non-trivial local stable manifolds for the Poincare map '■ Mc Mc- We will use a 
classical stable manifold theorem for hyperbolic compact set; so, we will be left to prove that the compact 
set C n Mc is a hyperbolic set for the Poincare map ^c- The key point is to understand the behavior 
of the orbits of the Wainwright-Hsu vector field close to a point p ^ C. Roughly speaking, we will 
prove the following: when an orbit of passes close to a point p G C, the distance of this orbit to the 
mixmaster attractor A = ICUBji decreases super-linearly, while the drift of this orbit tangentially to A 
is as small as wanted. We do not have any precise control of what happens to the orbits of Xq far from 
the Kasner circle, but we do not need to. Indeed, the duration of the excursion of the orbits of Xjg outside 
any given neighborhood of /C is universally bounded. A consequence is that everything that happens 
far from /C is dominated by the super-linear contraction of the distance to A that occurs when an orbit 
passes close to C. This will be enough to obtain the desired hyperbolicity result. As already explained, 
the reason why we need to consider an aperiodic subset of /C is purely technical: close to a point of /C 
which is not preperiodic under the Kasner map, we have a "nice" coordinate system which is convenient 
to study precisely the behavior of the orbits of X^. 

Let us now explain the organization of the paper. Consider a vector field X on some manifold and 
a point p such that X{p) = 0. We say that X satisfies Sternberg's condition at p if the non purely 
imaginary eigenvalues of DX{p), counted with multiphcities, are linearly independent over Q. Takens 
has proved a generalization of Sternberg's theorem which states that, if X satisfies Sternberg's condition 
at p, there exists a local "linearizing" coordinate system for X on a neighborhood of p. A precise 
statement of this theorem will be given in section 2. The purpose of section 3 is to apply Takens' theorem 
to the Wainwright-Hsu vector field. It will be proved that the Wainwright-Hsu vector field Xjs satisfies 
Sternberg's condition at some point p of the Kasner circle if and only if p is not preperiodic under the 
Kasner map. In section 4, the "linearizing" coordinate system provided by Takens' theorem will be used 
to study the behavior of the orbits of the Wainwright-Hsu vector field close to a point of the Kasner circle 
which is not periodic under the Kasner map. Now consider a closed forward-invariant aperiodic subset C 
of the Kasner circle /C is considered, and denote by C the union of C and all the type II orbits connecting 
two points of C in B'^. A Poincare section Mc adapted to C will be constructed in section 5; we will 
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denote by the corresponding Poincare map. The existence of non-trivial stable manifolds for the 
corresponding Poincare map '■ Mc Mq will be proved in section 6, using the results of section 4. 
The proof of theorem 1.3 will be completed in section 7. Finally, the last section of the paper is devoted 
to the proof of proposition 1.6. 

Acknowledgements. 1 would like to thank Lars Andersson for some stimulating discussion, and for his 
encouragements to write the present paper. 

2. TAKENS' linearization THEOREM 

Let X be a C°° vector field on some manifold A4, and p be a point in M. such that X{p) = 0. The 
linear space TpM admits a unique decomposition 

TpM = E'®E^eE'' 

where E^, E'^, E^ are Z)X(p) -invariant linear subspaces, the eigenvalues of DX{p)^Es have negative 
real parts, the eigenvalues of DX{p)\E': are purely imaginary, and the eigenvalues of -DX(j»)|£;u have 
positive positive. We denote by s, c and u the dimensions of the Unear subspaces E^, E'^ and E'^. 

Definition 2.1. The vector field X satisfies Sternberg 's condition at p if the eigenvalues of DX{p)\es^e" 
(counted with multipHcities) are linearly independent over Q. 

F. Takens as proved the following generaUzation of the classical Sternberg linearization theorem: 

Theorem 2.2 (see [10], page 144). Assume that X satisfies Sternberg's condition at p. Then, for every 
r > 0, one can find a neighborhood U and a C" coordinate system {xi, . . . ,Xs,yi, ■ ■ ■ ,yc, zi, . . . , Zy) 
on U centered at p, such that, in this coordinate system, X reads: 

^ d d ^ d 

(3) ^=Y^ ('iAyi^---^yc)xjQ^ + ^4>i{yi.---,yc)-Q-+Yl Kjiyu---^yc)zjQ^ 

i,j=l * 1=1 i,j=l ' 

where the eigenvalues of the matrix (aij(0, . . . , 0))i<ij<s have negative real parts, the eigenval- 
ues of the matrix (f^(0, • • • , 0))^^. .^^ are purely imaginary, and the eigenvalues of the matrix 
{bij{0, . . . , 0))i<ij<„ have positive real parts. 

Note that, in general, the size of the neighborhood U does depend on the integer r, and it is not 
possible to find any local coordinate system centered at p such that (3) holds. In this sense, Takens' 
theorem is not a true generalization of Sternberg's theorem. 

Also observe that the name "Takens' linearization theorem" is sUghtly incorrect: indeed, the vector 
field X is not hnear in the (xi, . . . , x^, yi, . . . , yc, zi, . . . , Zu) coordinate system. Nevertheless, X de- 
pends linearly on the coordinates xi, . . . ,Xs and zi, . . . ,Zu- Also note that the submanifold defined by 
the equation (yi, . . . ,yc) = (0, . . . , 0) is invariant under the flow of X, and that the restriction of X 
to this submanifold is linear. For {(^i,. . . ,(c) 7^ (0, • • • , 0), the submanifold defined by the equation 
(yi) ■ ■ • ) Z/c) = (Cii • ■ ■ 5 Cc) is not invariant under the flow of X, but the projection of X on the tangent 
space of this submanifold is linear. 

Of course, equality (3) together with the signs of the real parts of the eigenvalues of the matrices 
(aij(0, . . . , 0))i<ij<s, (g^(0, . . . , 0))i<ij<c (6ij(0, . . . , 0))i<ij<„ implies that: 
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- the vectors ^{p), ■ ■ ■ , span the Unear subspace E^, 

- the vectors -^{p) , . . . , {p) span the linear subspace E*^, 

- the vectors ^(p); ■ • • > ^(p) span the Unear subspace E'". 

3. Linearization of Wainwright-Hsu vector field near a point of the Kasner circle 
which is not preperiodic under the kasner map 

The purpose of this section is to apply Takens' theorem to the Wainwright-Hsu vector field at a point 
of the Kasner circle. For this purpose, we will need to characterize the points p on the Kasner circle 
such that the Wainwirght-Hsu vector field Xq satisfies Sternberg's condition at p. We will see that these 
are exactly the points p e K which are not preperiodic under the Kasner map /. In order to relate the 
arithmetic properties of the eigenvalues of the derivative of DXq{p) and the behavior of the orbit of p 
under /, we will use the so-called Kasner parameter. 

3.1. Kasner parameter. Let q = (Si, S2, S3, 0, 0, 0) be a point of the Kasner circle. The Kasner 
parameter of q is the unique real number u = u{q) G [1, 00] which satisfies the following equaUty: 

(4) (2 + Si) (2 + S2) (2 + S3) = I J ■ 

The map q 1— u{q) is not one-to-one. Nevertheless, the point q is characterized by its Kasner parameter 
up to permutations of the coordinates Si, S2, S3. More precisely, equahty (4) together with the equation 
of the Kasner circle (2) imply that: 

(5) (2 + E.,2 + 2E..2 + 2E3) ""J", (r^.j^.^^). 
^ ^ permutation Xl + u + W^ l + u + W^ l + u + W^J 

Note that n = 00 if and only if q is one of the three special points Ti, r3. The main advantage of the 
Kasner parameter is the fact that the Kasner map / admits a nice expression in terms of this parameter: 
for every g G /C, one has u{f{q)) = f {u{q)) where / : [1, +00] — > [1, 00] is defined by 

00 if II = 1 or 00 



(6) fiu) = { 



u-1 ifu>2 



ifl<u<2 



(see, for example, [3]). 



3.2. Characterization of the points of the Kasner circle where Sternberg's condition is satisfied. 

The proposition below gives a necessary and sufficient condition for the Wainwright-Hsu vector field Xq 
to satisfy Sternberg's condition at a point p G /C, in terms of the behavior of the orbit of p under the 
Kasner map. The hypothesis of our main theorem 1.3 comes directly from this condition. 

Proposition 3.1. Let p be a point of the Kasner circle K which is not one of the three special 
points Ti, T2, Ts. The three following conditions are equivalent : 

(1) the vector field Xq satisfies Sternberg 's condition at p ; 

(2) the Kasner parameter u{p) is neither a rational number, nor a quadratic irrational number ; 
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(3) the orbit ofp under the Kasner map f is not preperiodic. 



Proof. Denote by (Ei, E2, S3, 0, 0, 0) the coordinates ofp. Since pis not one of the three special points, 
the derivative DX^{p) has two distinct negative eigenvalues, one zero eigenvalue, and one positive 
eigenvalue. The three non-zero eigenvalues of Z)Xg(p) are equal to — (2+Si), — (2+S2) and — (2+S3). 

Let us prove the equivalence between (1) and (2). The vector field satisfies Sternberg's condition 
at p if and only if the real numbers —(2 + Ei), —(2 + E2) and —(2 + E3) are linearly independent 
over Q. Using formula (5), one sees that this is equivalent to the fact that the real numbers u{p), 1 + u{p) 
and u{p){\ + u{p)) are independent over Q. Clearly, this is equivalent to the fact the real number u{p) 
is neither a rational number, nor a quadratic irrational number. 

Now, let us prove the equivalence between (2) and (3). Recall that, for every q on the Kasner circle, 
one has u{f{q)) = f{u{q)) where / : [1, +00] [1, oo] is given by (6). Observe that both the set of 
rational numbers and the set of irrational numbers are invariant under /. So we can treat the case where 
u{p) is rational and the case where u{p) is irrational separately. 

First consider the case where u{p) is rational. Then it is very easy to prove that the orbit of u{p) under 
/ "ends up" at 00. Now recall that u{q) = 00 if and only if q is one of the three special points. Hence the 
orbit of p under / "ends up" at one of the three special points. In particular, the orbit of p is preperiodic. 

Now consider the case where u{p) is irrational. Looking again at (6), one sees that the orbit of u{p) 
under / returns an infinite number of times in the interval (1, 2]. Let k : (1, 00) \ Q — > N \ {0} be the 
return time function of / in the interval (1,2], and F : (1,2]\Q— >■ (l,2]\Qbe the first return map of 
/ in the interval (1,2], that is 

k(u) = inf{n > such that /"(u) G (1, 2]} and F{u) = /'=(")(t 

Then, 



k{u) 



1 and F{u) = <^ \ + 1, 



u 



where [xj is the integer part of x and {x} := x — [xj is the fractional part of x. The point p is 
preperiodic under the Kasner map / if and only if either u{p) is a preperiodic under /, that is if and only 
if jK'^ip)) {u(jp)) is preperiodic under F. Now observe that F is just the Gauss map u 1—^ j;^} conjugated 
by the translation u+1, and that /c(m — 1) + 1 is the first term of the continuous development fraction 
of u. On the one hand, the preperiodic points of the Gauss map n { ^ } are exactly the real numbers 
u G (0, 1] such that the sequence of integers which appear in the continuous fraction development of u 
is preperiodic. On the other hand, it is well-known that the continuous fraction development of u € M is 
preperiodic if and only if ti is a quadratic irrational number. This shows that the orbit of p is a preperiodic 
under the Kasner map / if and only if u{p) is a quadratic irrational number. □ 



3.3. Linearization of the Wainwright-Hsu vector field. According to proposition 3.1, if p is not a 
preperiodic point for the Kasner map, the hypotheses of Takens linearization theorem 2.2 are satisfied by 
at p. This theorem provides us with a local coordinate system on a neighborhood of p in which Xq 
is "almost linear": 

Proposition 3.2. Let pbe a point of the Kasner circle K which is not preperiodic under the Kasner map 
( in particular, p is not a special point). Then there exists a neighborhood U ofp in B and a coordinate 
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system (xi, X2, y, z) on U, centered at p, and such that, in this coordinate system, Xq reads: 

d d d 

(7) XB{xi,X2,y,z) = \\{y)xi— + \l{y)x2^^ + X'{y)z— 

where Af (y) < A^y) < and A"(y) > 0/or every y. 

Remarks 3.3. Let [/ be a neighborhood of the point p and (xi, X2, y, z) be a coordinate system on U 
centered dXp, such that Xq satisfies (7) with \\{y) < X^^iv) < A"(y) > 0. Then: 

(i) The vector field Xjs vanishes on the one-dimensional submanifold {xi = X2 = 2; = 0} and 
nowhere else. It follows that this submanifold is the intersection of the Kasner circle with U : 

/Cn [/ = {xi = X2 = z = 0}. 

(ii) For every y, the real numbers Af (y), A^y), A"(i/) are the three non-zero eigenvalues of the 
derivative of Xb at the point of coordinates (0,0, y,0) G /C. Recall that this derivative also 
has one zero eigenvalue (corresponding to the direction of the Kasner circle). 

(iii) For every (, the three-dimensional sub-manifold {y = ^} is invariant under the flow of Xq and 
the restriction of Xq to this submanifold is linear. 

(iv) The three-dimensional submanifolds {xi = 0}, {x2 = 0} and {z = 0} are invariant under the 
flow of Xi^, and contain K, r\U. It follows that these submanifolds coincide up to permutation 
with the submanifolds {A^^i = 0}, {A'^2 = 0} and {N^ = 0}. As a consequence, the two- 
dimensional submanifolds {xi = X2 = 0}, {xi = 2; = 0} and {x2 = 2; = 0} coincide up to 
permutation witii the submanifolds {A^i = N2 = 0}, {A^i = A^3 = 0} and {N2 = N3 = 0}. In 
particular, 

{Bii n /C) n = {xi = X2 = 0} u {xi = z = 0} u {x2 = 2 = o}. 

(v) The right-hand side of (7) is unchanged if one replaces xi (resp. X2 and z) by — xi (resp. by 
— X2 and —z). Therefore, we may assume that 

^+ n [/ = {xi > , X2 > , z > 0}. 

(vi) On the one hand, according to item (i), the metric induced on the one-dimensional submanifold 
/C n [/ by the riemannian metric (dxi)^ + {dx2)^ + (dy)^ + (d^)^ is simply (dy)^. On the other 
hand, the right-hand side of (7) is unchanged if one replaces y by ip{y) where 93 : R ^ R is 
any diffeomorphism such that ip{0) = 0. Therefore, up to replacing the coordinate y by (p{y) for 
some appropriate diffeomorphism (p, we may assume that the metrics induced on the piece of 
circle ICnU hy the riemannian metrics (dxi)^ + ((ix2)^ + (dy)"^ + (dz)"^ andby the riemannian 
metric h = ((iSi)^ + {d^2)^ + {dll^f + {dN^f + {dN2f + {dN3)2 coincide. 

Proof of proposition 3.2. The derivative DXb{p) has two negative, one zero, and one positive eigen- 
value. According to proposition 3.1, the vector field X^ satisfies Sternberg's condition at p. There- 
fore, a crude application of Takens' theorem 2.2 implies that there exists a local coordinate system 
(xi, X2, y, z) on a neighborhood [/ of p in B, centered at p, such that : 




for some real valued functions 4>, ai,i, 01,2, 02,1, «2,2) i> defined on a neighborhood of in R. Moreover, 
the eigenvalues of the matrix (aij(O)) are negative, and 6(0) is positive. Replacing i7 by a smaller 
neighborhood of p if necessary, we can assume that the three Taub points Ti,T2, Ts are not in U. 
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Now, (8) implies that the curve {xi = X2 = = 0} is the only curve in U containing the point p, 
invariant under the flow of Xb, and such that DXb{p) vanishes on the tangent space at p of this curve. 
Hence, the curve {xi = ,T2 = ^ = 0} has to be the intersection of the Kasner circle /C with U. Since Xb 

vanishes on /C, it follows that = 0. 

For C, small enough, let be the point of coordinates (xi, X2, y, z) = (0, 0, C, 0). This is a point of the 
Kasner circle /C, which is not a Taub point. Hence, the derivative DX^qq) has four distinct eigenvalues : 
two distinctive negative eigenvalues Af {Q < A2 (C) < 0, one zero eigenvalue associated to the direction 
of the Kasner circle, and one positive eigenvalue A"(C) > 0. The set (y = C) is a three-dimensional 
manifold, transversal to the Kasner circle. Looking at (8), we see that this three-dimensional submanifold 
is invariant under the flow of Xb, and that the restriction of Xb to this invariant manifold is linear. This 
shows that 6(C) = A"(C), and that Af(C) and A2(C) are the eigenvalues of the matrix (at,j(C))ij=i,2- 
Since Af {Q and A2(C) are distinct, there exists a hnear change of coordinates (xi, X2, z) {x\,X2,z) 
on the submanifold {y = Q, so that 

Xe(xi,x2,C,^) = Af(C)xi^ + Ai(C)x2^ + A"(C)^^. 

Since eigenvalues and eigendirections of the point depend in a smooth way of C,, one may perform the 
above change of coordinates simultaneously for every C„ and get a coordinate system (xi, ^2, y, z) 
defined on U, such that 

XB{xi,X2,y,z) = \{{y)xi-^ + \l{y)x2^ + y{y)z-^. 

0x1 0x2 oz 

The proposition is proven. □ 



3.4. Choice of a linearizing coordinate system. From now on, for every point p of the Kasner cir- 
cle which is not preperiodic under the Kasner map, we fix a neighborhood of p in B, and a 
local coordinate system (xf, y^, z'p) on V^, centered at p, such that, in this coordinate system, the 
Wainwright-Hsu vector field Xb reads: 

(9) Xb {xl,xly,zn = A?(/)x?^ + Ai(/)x^^ + X^iy^zP^ 

with Af (yP) < A|(y^) < and A"(y^) > 0. According to item (ii) of remarks 3.3, the real numbers 
Xfiy^), A|(y^), A"(y^) are the non-zero eigenvalues of the point (0,0,1/^,0). According to items (i), 
(iv), (v) of remarks 3.3, up to changing the sign of the coordinates x^, x?^ and z^, we may (and we will) 
assume that 

(10) /Cnt/*' = {xf = xf = z = 0} 

(11) (HnU/C)n[/f = {x? = xf = 0}u{x? = zP = 0}u{xf = zP = o}, 

(12) n = {x? > , xf > , > 0} 

Note that the derivative DXB{q) has three non-zero eigenvalues for every point q e KCi U^; this shows 
that the neighborhood is disjoint from the three special points. We will consider the riemannian 
metric on defined by 

(13) / := (dx?)' + (dxf )' + {dyPf + {dz^) 2. 
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According to item (vi) of remarks 3.3, up to replacing the coordinate by <y?(j/^) for some appropriate 
diffeomorphism (p, we may (and we will) assume that induces the same metric on the piece of Kasner 
circle /C n i/f as the riemannian metric h = {dEif + (dEa)^ + {d^sf + {dNi)^ + (ciiVs)^ + {dN^f. 



4. DULAC MAP FOR WAINWRIGHT-HSU VECTOR FIELD NEAR A POINT OF THE KASNER CIRCLE 
WHICH IS NOT PREPERIODIC UNDER THE KASNER MAP 

Let p be a point of the Kasner circle which is not preperiodic for the Kasner map. The purpose of the 
present section is to analyze the behavior of the orbits of the Wainwright-Hsu vector field Xq close to p. 
More precisely, we want to consider an orbit of which passes close to p, and to study the evolution 
of the distance from this orbit to the mixmaster attractor A = ICU as well as the drift of this orbit in 
the direction tangent to the mixmaster attractor. 

4.1. The flow of Xq inside U^. We consider the neighborhood U^, and the coordinate system 
{xi,x^,yP, z'P) defined in subsection 3.4. Using the expression (9), one can calculate explicitly the 
time tmap X^ of the flow of the Wainwright-Hsu vector field Xb in the (xf , X2,y^, z^) coordinate sys- 
tem. It reads 

(14) Xl (x?, , yP, zP) = (x\e^'^^y'^' , e^^^^'')* , / , e^"^^'')*) . 

Of course, this expression is only valid as long as the orbit remains in the neighborhood U^. 

4.2. The box V^. Now, we fix some constants a < j5 and 7 > 0, and we consider the subset 
of [7*^1^+ defined by: 

(15) = yP(a, /3, 7) = {0 < X? < 7 , < < 7 , < z?' < 7 and a < < 

We assume that a, /3, 7 are small enough, so that is contained in the interior of U^. 

Remark 4.1. The set is a neighborhood of the point p in 5+ if and only if a < < /3. It is important 
to note that the results of the present section are valid even if p is not in V^. 

In the {x\, y^, z^) coordinate system, the set is a 4-dimensional box, i.e. the cartesian product 
of four closed intervals. The boundary of is made of eight faces. Three of these eight faces will play 
an important role in the remainder of the paper: 

(16) Mf := n {x? = 7} Mf := n {xf = 7} := f\ {z^ = -f} . 

Looking at (9), we notice that Xq is transversal to Mf , and N^, and is tangent to the five other faces 
of V^. Moreover, we notice that Xb is pointing inward along Mf and Ml"; it is pointing outward 
along A^^*. It follows that: 

• an orbit of X^ can enter in by crossing either the face Mf or by crossing the face M^* ; 

• an orbit of Xs can only exit by crossing the face A'^*'. 
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4.3. Behavior of type 11 orbits. We will study the behavior of the orbits of Xjs in V^. First, we focus 
our attention on type n orbits. We want to understand which type 11 orbit intersect the hypersurfaces 

(with boundary and comers) Mf, Ml* and N^. Recall that every type II orbit O of Xq is a heteroclinic 
orbit connecting a point a{0) G /C \ {Ts, T2, T3} to the point uj{0) = /(a(C')) G /C \ {T2, Ts, T3}. 

Proposition 4.2. Let q be a point B^, and denote by O the orbit of q. Let a{q) = a{0) denote the 
unique a-limit point ofO, and oj{q) = uj{0) = f{a{0)) denote the unique oj-limit point ofO. 

(1) The orbit O intersects the hypersurface if and only if the point a{q) is in V'^. 

(2) The orbit of O intersects the hypersurface U Ml if and only if the point Lo{q) is in V^. 

Proof. We prove the first statement; the second one follows from similar arguments. Of course, we will 
work in the (xf , X2, y^, z^) coordinate system. According to (10), (11), (15) and (16), 

ICnVP = {,T? = a;^ = 0, a</</3, zP = 0}, 
BiitlNP = {x? = = 0, a < y^' < /3, = 7}. 

Suppose that the orbit O intersects the hypersurface at some point q = (0, 0, C, 7), with a < ( < 13. 
Then, according to (14), the past orbit of q is contained in V^, and converges to the point a{q) = 
(0, 0, C, 0) G /C n VP. In particular, the point a{q) is in V^. Conversely, suppose that the point a{q) is 
in VP. Then a{q) = (0, 0, 0) for some ( G [a, (3]. Using again (14), we see that the only orbit of 
in VP converging towards the point (0, 0, (, 0) as t — —00 is the curve t (0, 0, (, e^"^'')*). Hence, the 
orbit O intersects the hypersurface A^^ at the point ^ = (0, 0, C, 7). □ 

This proposition allows to define two maps 

aP: B+nNP ^ }Cn VP ojP : n (Mf U Mf ) ICnVP 

q I — > a{q) q 1 — > uj{q) 

The map is one-to-one (there is only one type II orbit in which "starts" at a given point of 
/C \ {T2, Ts}), whereas the map ujP is two-to-one (there are two type II orbits in B^ which "arrive" at 
a given point of /C \ {r2, T2, T3}). The restriction of ujP to B^^ D Mf (resp. to 6+ n M|) is one-to-one. 

Proposition 4.3. The maps aP and ujP are local isometrics with respect to the metrics induced on 
BfjHNP, B^jniMfuM^) and ICnVP by the riemannian metric gP = {dx^)'^ + {dxlf+{dyPf + {dzP)2. 

Proof. The proof of proposition 4.2 shows that, in the (x^ x^, yP, zP) coordinate system, the map aP 
reads aP{<d, 0, yP, 7) = (0, 0, yP , 0). Similarly, the map uP reads a;P(7, 0, yP , 0) = (0, 0, yP, 0). □ 

4.4. The Dulac map : Mf NP. Now, we want to study the behavior of arbitrary orbits of Xb 

which enter in VP. Let g be a point on the face Mf . Denote by (7, X2, yP, zP) the coordinates of q. If the 
> (which is typically the case if g' G 6ix)» then (14) shows that the forward orbit of q will eventually 
exit VP by crossing the face NP. If = (which is typically the case if g G B-^^), then (14) shows 
that the forward orbit of q will remain in VP forever. It will converge towards the point uj{q) ^lCf\ VP. 
According to proposition 4.2, the heteroclinic orbit ©(^(q) will eventually exit VP, by crossing the 

face NP. So, we may define a map : Mf — > NP as follows : 

• if > then A^(g) is the first intersection point of the orbit of q with the hypersurface AT*' ; 



14 FRANCOIS BEGUIN 

• if z^' = then A^(^) is the first (and unique) intersection point of the type II heteroclinic orbit 
Ouj{q),f{oj{q)) with the hypersurface NP. 

We call a Dulac map since it is the exact analog, in our situation, of the classical Dulac maps used to 
study planar vector fields. Formula (14) show that, in the (xf , x^, y^, z^) coordinate system, the map 
reads : 




(17) A{{j,x'^,yP,zn = |7-{-) '^2-{-) ,/,7|if^^>0 

(18) Af(7,xf,y^0) = 

Remark 4.4. Given a point q in Mf such that z'^ (q) > 0, one can consider the exit time of q, that is the 
real number t{q) such that A^(g) = X^^'i\q). Using (14), it is easy to see 



i(7,x^/,.^) = -^log(^) 



For every q G Mf, we decompose TqM^ as a direct sum of two linear subspaces © where 

(19) := ^-^W^^-^W 

(20) := M.A(,). 

Similarly, for every r G Np, we decompose TrN^ as a direct sum of two linear subspaces © 
where 

(21) := K.A(,)eM.^(r) 

(22) := M. — (r). 

We can now state the properties of the Dulac map A^ which will be the core of our proof of theorem 1.3: 

Proposition 4.5. The Dulac map Af : Mf NP is C^. Moreover, for every point q G Mf n {zP = 0}, 
the derivative of the map A\ at q satisfies: 

• DA\{q).v = for every vector v G F^; 

• DA^(q) maps F^ on G^p^^y and \\DA^{q).v\\gP = \\v\\gpfor every vector v G F^. 

This proposition roughly says the following: when an orbit of passes close to the point p G /C, the 
distance from this orbit to the mixmaster attractor A = ICU Bn is contracted super-linearly (this distance 
is measured by the coordinates x^, x^ and zP), whereas there is no drift in the direction tangent to the 
attractor (this drift is measured by the coordinate yP). The key point of the proof of proposition 4.5 is the 
following elementary observation: 



Lemma 4.6. For every yP, we have 



\^{yP) 



> 1 and 



A"(yP) 



> 1. 
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This lemma says that, at every point g of /C \ {T2, T2, T^}, the positive eigenvalue of the derivative 
DX{q) is dominated by the contracting eigenvalues. 

Proof of lemma 4.6. Fix y^, and denote by q be the point of coordinates (0, 0, y^, 0) in the {x\, X2, y^, z^) 
coordinate system. Denote by (Si , S2 , S3 , 0, 0, 0) the coordinates of q in the Wainwright-Hsu coordinate 
system, and by u the Kasner parameter of q. The real numbers Af (y^), X2{y^), A"(yP) are the three non- 
zero eigenvalues the derivative DXjs{q)- Hence, these numbers are equal up to permutation to — (2+ Si), 
-(2 + S2), -(2 + S3). Using (5) and the inequalities \\{yP) < A|(yP) < < A"(yP), we deduce that 

(ASfe-) , K(yn , A"(/)) = (=5^^ , =^^, , ^ , 

The lemma follows since > 1 and (1 + u) > 1 for every n > 1. □ 

Proof of proposition 4.5. The fact that is (and even analytical) in restriction to n {z^ > 0} 
is an immediate consequence of formula (17). The fact that Aj* is on n {z^ = 0} follows 
from(17), (18) and lemma 4.6. The same ingredients actually show that, for every point g G MPr]{zP = 
0}, 

= = andi^A?(g). A(g) = |^(A?(g)). 

The proposition follows. □ 
Remark 4.7. We do not know if A^ is C^~^^ for any given £ > 0, unless we have some a priori lower 
bounds for the distance between the ratios x^{yp) XHv^) ^' reason why, in the state- 

ment of theorem 1.3, we cannot guarantee that the stable manifold W^{q) contains a C^+^-embedded 
disc for any given s. Actually such an e exists for every q but it does depend on q, and tends to if g 
approaches one of the three special points tends to 0. 

4.5. The Dulac map Aj : — > N'^. The coordinates and X2 play similar roles in the expression of 
the vector field Xq and in the definition of V^. So, we have to consider a second Dulac map Ag : — >■ 
NP defined as follows. Let g be a point on the face Mf . Denote by (7, X2, y^, z^) the coordinates of q. 

• if z^{q) > then A2(g) is the first intersection point of the orbit of q with the hypersurfaceA'^^ ; 

• if zP{q) = then A2{q) is the first (and unique) intersection point of the type II heteroclinic orbit 
Ooj{q),f{uj{q)) with the hypersurface NP. 

In the (xf , x^,y''',zP) coordinate system, the Dulac map A2 reads: 

(24) A^(x?,7,y^0) = (0,0,y^7) 

For every q G Mf , we will write T^Mf as a direct sum of two linear subspaces © F^ where 



(25) := M.^(g)©M. — (g) 



(26) F^ := M.A(g). 
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Then, we can summarize the key properties of the map A2 as follow : 

Proposition 4.8. The map : Mf -> A^^ is C^. Moreover, for every g G M| n {z^* = 0}, the 
derivative of of the map A2 at q satisfies : 

• DA2{q)-v = for every vector v G ; 

• DA2(q) maps on G^p^^y and \\DA2{q).v\\gP = \\v\\gp for every vector v in F^. 

5. Construction of a "Poincare map" associated to a closed forward-invariant 

APERIODIC SET OF THE KASNER CIRCLE 

From now on until the end of section 7, we consider a closed forward-invariant aperiodic subset C of 
the Kasner circle IC. 

Observe that C is necessarily totally discontinuous. Indeed the points of K, that are preperiodic under 
the Kasner map / are dense in /C (this can be proved in several different ways; this follows for example 
the equivalence 2 44> 3 of proposition 3.1). We shall denote by C the union of C and all the type n orbits 
connecting two points of C: 

We want to prove that, for every point q & C, the stable manifold W'^{q) contains a 3-dimensional disc 
(see definition 1.2 and theorem 1.3). To this end, we will consider a kind of "Poincare section" for C, 
and study the Poincare map associated with this section^ 

For every jo € C, we consider a neighborhood of p, and a local coordinate system {x^, y^, z^) 
on U'P, as in the previous section. For each p ^ C,^e choose Qq < < /3q and 7o > small enough, 
so that 

:= V^ial^l^l) = {0 < X? < 7o^ < .T^ < 7^ < < 70" and al<y^< 

is contained in the interior of U^. Up to slightly modifying ckq, /3q , 7q we can assume that the boundary of 
Vq is disjoint from C (i.e that the points of coordinates (0, 0, a, 0) and (0, 0, /3, 0) in the (xf , x^, y^, z^) 
coordinate system are not in C): this is possible since C is totally discontinuous. Observe that Vq* is a 
neighborhood of p, since < < /3q (see remark 4.1). 

Since C is compact, one can find a finite number of points pi , . . . , p„ € C such that the neighborhoods 
Vq*^ , • • ■ , Vq" cover C. Now, we modify these neighborhoods in order to make them pairwise disjoint: 

• we set (aPi , /3pi , 7^1 ) := ("o' ' /^o' > 7o' )' ^nd V^^ := {uP^ , j3P^ , 7^^ = V^^ ; 

• then, we can find some constants a2,/32)72 such that V^'^ := VP'^{a2j 152,12) is contained in 
yP2 ^^ ypi ^ ^jj^ gypjj ^jj^^ (J p ^yP2 ^ ypi ) contained in the interior of • 

• then, we can find some constants q;3,/33,73 such that V^^ := T^^'^(a3, /Ss, 73) is contained in 

\ {VP^ U VP^), and such that C n {V^^ \ {VP^ U VP^)) is contained in the interior of VP^; 

• etc. 

^The set C cannot admit a true Poincare section, since it contains some singularities of Xb (namely, the points of C). 
Nevertheless, we wiU consider a hypersurface N such that every type II orbit in C intersects N transversally. The hypersurface 
N will play the role of a Poincare section. 
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At the end of this process, we get n pairwise disjoint domains V'^^ , . . . , V^", such that C is contained in 
the interior of V^'^ U ■ ■ ■ U F*'". For each i, V^^ is contained in the interior of UP\ and there are some 

constants ctj, /3j, 7j such that V^^ = V^^{ai, jii, 7j). Hence, the result of section 4 apply to It may 
happen that, for some i, the point is not in V^^ {i.e. that or fii is non-positive), but we do not care. 

Now, we denote 





:= V^^ 


U • 


• u yp", 


Mf 


:= Mf 


U • 


• • U Mf " 


Mf 


:= Mf 


u • 


• • U " 




:= Mf 






:= iV^i 


U • 


• • U iV^". 



Then 1/*^ is a neighborhood of C in . The hypersurfaces M{^, M2 and A^*^ are transverse to X^- An 
orbit of Xq can only enter in V'~^ by crossing M*^ = U Mg^, and can only exit by crossing A'^'-^. 
Moreover, according to proposition 4.2, we have the following important properties : 

Proposition 5.1. Every type II orbit whose uj-Umit point is in C intersects M^. Every type II orbit whose 
a-limit point is in C intersects N^. 

We will see M*^ as a kind of "Poincare section" for C. Let us define the "Poincare map" $ associated 
to this section. First, we consider the "Dulac map" 

defined by Aj^p^ = Af and Aj^^p^ = Af. For q e M^: 

• if the forward orbit of q exits V*^ by crossing N*^ (which is typically the case if g G Bix), then 
A*-^(g) is by definition the first intersection point of the orbit of q with the hypersurface ; 

• if the forward orbit of q remains in forever (which is typically the case for every q G 
then A'^{q) is the first intersection point of the type 11 heteroclinic orbit Cu)(q),/(t^(5)) with the 
hypersurface A^*^. 

Now, we consider the "transition map" 

partially defined as follows. Given a point q in N"^ , if the forward orbit of q re-enters in , then {q) 
is the first point of this forward orbit of q which is in V'-' (this point is automatically on the hypersurface 
M*^); otherwise 6*" (g) is not defined. The "Poincare map" associated with the section is by 
definition the product of the "Dulac map" A^ and the "transition map" 

^C' qC oA^ •.M'^ ^ M'^. 

In the next section, we will study the dynamics of the Poincare map 0*^. For this purpose, we will use 
a riemannian metric on on B such that, for i = 1, . . . , n. 
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6. Stable manifolds for the Poincare map associated to a closed 

FORWARD-INVARIANT APERIODIC SET OF THE KASNER CIRCLE 

The purpose of this section is to prove that, for every point q & C Ci , the stable manifold of q for 
the "Poincare map" contains a two-dimensional disc. To this end, we will prove that C n is a 
hyperboUc set for the map and we will use a classical result on stable manifolds for hyperboUc sets. 

Definition 6.1. Let (M, g) be a riemannian manifold and $ : M — > M be a map. A hyperbolic 
set for the map $ is a compact ^-invariant subset C of M such that, for every q E C, there is splitting 
TqM = ® which depends continuously on q and such that, for some constant ji G (0, 1) and 
nu> 1: 

(27) D^{q).F^ C and ^^(x).^!! < ii\\v\\ for every q^Czndv eF^ 

(28) D^{q).F^ = and ||£>$(g).u|| > v\\v\\ for every g G C and u e F^. 

The dimension of the vector space F^ is called the index of C. The constant /x is called a contraction 
rate of $ on C. 

Theorem 6.2. (see e.g. [6, page 167] j Let ^ : M ^ M be a map of a manifold M, and C be a 
compact subset of M which is a hyperbolic of index sfor the map Then, for every e small enough, for 
every q E C, the set 

Wl{^, q):={r(^M\ J/5?($"(r), $"(r)) < e for every n > 0} 

is a embedded s-dimensional disc, tangent to F^ at q, depending continuously on q (for the 
topology on the space of embeddings). Moreover, if p, is a contraction constant for # on C, then there 
exists a constant k such that, for every e small enough, for every q E C and every r G q). 

We want to apply this theorem to the Poincare map : — > M'^. So we need to prove that 
d n M*^ is a hyperbolic set for Recall that M*^ = Mf U where Mf = Mf ' U • • • U Mf " and 
= U • • • U Mf". For every q G W^, we have already defined a splitting T^M'^ = F;^ ® F^ 
in section 4 (recall that M"^ = (Mf' U • • • U Mf") U (Mf^ U • • • U M|") and observe that q is not in 
(Mf 1 U • • • U Mf " ) n (Mf 1 U • • • U Mf " )). It remains to prove that satisfies (27) and (28) with respect 
to these splitting. For this purpose, we will use the decomposition of as a product : 

= e*^ o A*^. 

The behavior of the derivative of "Dulac map" was already studied in section 4 ; more precisely, we 
can rephrase propositions 4.5 and 4.8 as follows : 

Proposition 6.3. The map ISp : — > N'^ is C^. Moreover, for every q E C D M'^, the derivative of 

DA'^{q) : TgM T^c^q)N of the map ^ at q satisfies : 

• DA^ {q).v = ^ for every vector v G ; 

• DA^(^q) maps Fg on G^c(q). o,nd \\DA^ {q).v\\g^ = HuUgp/or every vector v in Fg. 

It remains to study the behavior of the derivative of the "transition map" Q*-" : — )■ . We recall 
that 0*^(g) is well-defined only if the forward orbit of q intersects M*-^. So our first task is to show that 
0^ is well-defined at least on a neighborhood of C n N'^ in N'^ . 
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Proposition 6.4. There exists a neighborhood V of C H in , such that, for every g G V, the 
orbit of q intersects after some time t{q) which depends in a way on q. The map 0*^ is well- 
defined and on V. Moreover, there exists u > 1 such that, for every q E C iV^, the derivative 
De^{q) : TqN^ TQC^q^M^ satisfies 

• L>e^(g).G« = F^c^q^ and \\DQC{q).v\\g > u\\v\\g for every v E 

Proof Consider a point g G C fl . By proposition 5.1, the orbit of q intersects at some point 
r G C* n M<^. Now, recall that : 

• N, , M2 are hypersurfaces with boundary that are transversal to the orbits of ; 

• V was chosen so that C is contained in the interior of V. This implies that C does not intersect 
neither the boundary of the hypersurface A^, nor the boundary of hypersurface M^' and It 
follows that C does not intersect n Mf^. Hence, q is in the interior of N, and r is in the 
interior of Mf or M^. 

This implies the existence of a neighborhood Vq of qinN such that, for every q' G V, the forward orbit 
of q' intersects M*^ after some time t{q') which depends in a way on q'. By definition of 9, for every 

q' G Vq, we have Q{q') := In particular, is well-defined and on Vq. This proves the 

two first statements of the proposition. 

Since Bn is invariant under the flow of Xg, the map Q'-^ maps Bn D N'-^ on Bn D M'-\ Now, observe 
that, for every q G Bn n N'-^, the direction is nothing but the tangent space of Bn n A^"-^ at q, and the 
direction -^0C(g) is nothing but the tangent space of Hjj fl at @'-^{q). This shows that d&'-^{q) maps 
on F^c(g) for every q e C H N'=' . 

We are left to prove the existence of a constant u > 1 such that ||Z)G*^(g).f ||g > for every 

q E C n N^' and every v G Gq. For this purpose, we will use the maps 

a : % n iV^ — ^ /C n and a; : % n ^ /C n . 

We recall that a maps a point r G Bn H A^*^ to the a-limit point of the orbit of r, and that to maps a 
point s G n M*-^ to the w-limit point of the orbit of s (see section 4). We also recall that a is a 
local isometry for the metrics induced by g on Bji D and /C fl F*^, and that w is a local isometry 
for the metrics induced by g on Hjj fl and /C fl V'~' (proposition 4.3). Finally, we observe that, for 

r G % n N^, 

a;(G^(r))=a;(r) = /(a(r)). 

The first equality is due to the fact that 0*^ (r ) and r are on the same orbit ; the second one is an immediate 
consequence of the definition of the Kasner map /. This shows that the last statement of proposition 6.4 
is equivalent to the following statement about the Kasner map : there exists a constant v > 1 such that, 
for every p EG and every v G TpK,, one has \\Df{p).v\\g > iy.\v\g. 

This last statement is an immediate consequence of the elementary properties of the Kasner map, and 
of our choice of the riemannian metric y*^. Indeed, the riemannian metric g^' was chosen so that it 
induces the same metric on /C n as the euclidean metric h = (c/Si)^ + (dS2)^ + {dT,^)'^ + {dNi)'^ + 
(diV2)^ + (dNs)"^ (see the end of subsection 3.4 and the end of section 5). And, as we already mentionned 
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in the introduction, since C is a compact subset of the Kasner circle /C which does not contain any of 
the three special points Ti, T3, there exists a constant i/'-^ > 1 such that, for every q e C and every 
V G Tp/C,wehave \\Df(p).v\\h > 2^*^111^11, where || • denotes the metric induced on /C by the euclidean 
metric h. □ 

Let U := $~^(V). Clearly, is a neighborhood of C n in M^. Combining propositions 6.3 
and 6.4, one immediately gets : 

Proposition 6.5. The Poincare map is well-defined and on U. The compact set C is a 
hyperbolic set for More precisely, there exists a constant v E {0,1) such that, for every q G 

• d^'~^{q).v = for every v G F^, 

• d^^{q).F^ = F^a^qy and \\{d^C {q))-\v\\g < u\\v\\g for every v e F^c^^y 

Proposition 6.5 shows that the map : M*^ — > M*^ and the set C Ci satisfy the hypotheses of 
the stable manifold theorem 6.2 (for any contraction rate ix > 0). This shows the existence of local stable 
manifold, with respect to the map <1>'^, for the points of C* H M*^ : 

Theorem 6.6. For every e small enough, for every q E C D , the set 

Wl{^^,q) ■- {r G I distg (($^)" (r), ($^)" (g)) < £ for every n> 0} 

is a -embedded disc of dimension 2 in , tangent to F^ at q, depending continuously on q in the 
topology. Moreover, for every constant n > 0, there exists another constant K such that, for every 
q e Ctl M^,for every r G W^{^'-', q) and every n > 

distg[{^'^Y{r),{^''Y{q))<K^x^. 



7. Stable manifolds for the Wainwright-Hsu vector field: proof of theorem 1.3 
We are left to prove that theorem 6.6 implies our main theorem 1.3. 

Consider a point q G C. Then heteroclinic orbit Oqj{^q) intersects the "Poincare section" at one 

and only one point, that we denote by q. Note that g G C fl M^. The set PF^ (<I>'-^, q) defined in the 
statement of theorem 6.6 is a C^-embedded two-dimensional disc in the three-dimensional hypersurface 
with boundary M'^ . This disc is tangent to F| at q. Since the two-dimensional submanifold ^vilo ^ B\i 
is not tangent to F| at q, this implies that q) fl ^Sjx contains a C^-embedded two-dimensional 

disc in . Moreover, this disc depends continuously on q. 

Proposition 7.1. For every point r in ($'-^, q) n Bjx-' 

(1) there is an increasing sequence of times (tn)n>o such that dist„{Xtp(r), f"^{q)) — > 0; 

(2) the Hausdorff distance between the piece of orbit {X^ir) ; tn <t < in+i} and the heteroclinic 
orbit Ofn(^qyfn+i(^q-^ tends to when n goes to -|-oo. 
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Proof. We first prove item 1. According to theorem 6.6, we have 

(29) distg f (cj.^)" (r), i^^Y {q)) 0. 



Together with the continuity of the flow of Xq, this shows the existence of a increasing sequence of real 
numbers {Tn)n>o such that 

(30) dist, (X- (($^)" (r)) ,u (($^)" (,-))) ^-^0. 

Since r G ;Bjx, there exists an increasing sequence of times (Q!n)„>o such that, for every n > 0, 

(31) ($^)"(r)=X^"(r). 
Since q G B^, we have, for every n > 0, 

(32) (($^)" {q)) = r{u{q)) = r+\q). 

For every n > 0, let tn+i ■= Tn + a„. Then (tn)„>o is an increasing sequence, and 

(33) dist,(4"+H0,r+'(9)) ^0. 
This completes the proof of item 1 . 

To prove item 2, we decompose the piece of orbit {X^{r) ; tn < t < tn+i} into three sub-pieces: 

• First, the piece of orbit going from X^{ro) to A*-^ ((<^>'-^)" (r)), contained in V. Formula (14), 
together with (33) shows that, for n large, this piece of orbit is close to the heteroclinic orbit 

^ particular, for n large, the point A*-^ (($'^)" (r)) is close to the heterocUnic 

orbit 0/n(g)Jn+l(g). 

• Then, a piece of orbit going from A*^ (($'^)" (r)) to ($*^)"^ (r), contained in ^ \ y. For n 
large, this piece of orbit is close to the heteroclinic orbit C'jn(g) jn+i(g). Indeed, for n large, 
the point A*-^ (($*^)" (ro)) is close to the heteroclinic orbit C'jn(q) jn+i(g), and if we write 

^ ^^(AC((*^)"(.))) (($C)- then t (A^ ((cD^)" (r))) depends con- 

tinuously on A*^ (($'-^)" (r)) (proposition 6.4) and thus is uniformly bounded. 

• Finally, a piece of orbit going from (^c-)""^' (r) to Xg"+'(r), contained in V. Formula (14) 
together with (33) show that, for n large, this piece of orbit is close to the heterochnic orbit 

This completes the proof of item 2 □ 
Corollary 7.2. 

w%q) = U U {^e (^"^^ {^''T (q)) n Bix) . 

t>On>0 

Proof. The inclusion of the set on the right hand side in W^{q) follows from proposition 7.1. The 
inclusion of W'^{q) in the set on the right hand side is an immediate consequence of the definition 

ofTy*(g). □ 



We can now complete the proof of our main theorem. 
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Proof of theorem 1.3. Fix > 0, and we set 

(34) D'^{q):= \J (W^/ ($^, n ^jx) • 

—ri<.t<.ri 

According to corollary 7.2, D^-^[q) is contained in Since W^{(^^ , q)) n is a -embedded 

2-dimensional disc in M*^, and since the orbits of Xq are transversal to M^, we get that -01^(5) is a 
-embedded 3-dimensional disc in Hjx- Since g depends continuously on q (proposition 4.3), and since 
W^{^'^ , q) depends continuously on q (theorem 6.6), the disc Dj^{q) depends continuously on q. □ 

Remark 7.3. The fact that (q) is an injectively immersed open disc which depends continuousluy 
on q (remark 1.4) almost follows from the same arguments. More precisely, theorem 6.6, corollary 7.2 
and the transversality of M*^ to the orbits of X show that W'^{q) is an increasing union of embed- 
ded closed discs which depend continuously on q. The only thing which remains to shows is that this 
increasing union of closed discs is an open disc; this is actually a consequence of the fact that the orbit 
of q under the Kasner map / is not periodic. 

8. Existence of closed forward-invariant aperiodic subsets of the Kasner circle: 

proof of proposition 1.6 

The purpose of this section is to prove proposition 1.6. This proposition should be quite obvious for 
people with some culture in dynamical systems. Indeed, the Kasner map / : /C — )■ /C is a degree —2 
map of the circle /C. This implies the existence of a continuous degree 1 map 77 : /C — > R/Z such that 
77 o / = m_2 o (j) where m : M/Z — > R/Z is defined by m,{6) = —29. Moreover, the Kasner map / 
is expansive : the norm of the derivative of / (calculated with respect to the metric induced on /C by 
the rimennian metric h) is strictly bigger than 1, except at the three special points Ti,T2, (where it is 
equal to 1). This implies that the map 77 : /C ^ M/Z defined above is one-to-one, that w / : /C — > /C is 
topologically conjugated to the map m : 9 ^ —26. Finally, it is well-known by experts that, for > 2, 
the union of all compact subsets of M/Z which are aperiodic for the map 9 ^ k9 is dense in M/Z. We 
now give a more detailed proof of the proposition for the readers who may not necessarily be familiar 
with low-dimensional dynamics. 

Proof. Recall that Ti,T2,T^ are the three Taub points on the Kasner circle K,. Let Ii,l2,h be the 
closures of the three connected components of /C \ {Ti, T2, Ts}, the notations being chosen so that Ti is 
not one end of Ii, T2 is not one end of I2, and T3 is not an end of I3. 

We consider the set S := {1, 2, 3}^ endowed with the product topology, and the shift map cr : S ^ S 
defined by a{ao, oi, 02, • • • )« = (ai, 02, 03, • • • ) (in other words, if a = (ai)igN € 5], then {a{a))i = 
Oj+i). Let So be the subset of E defined as follows : 

So = {a = (ai)igN G S such that aj+i 7^ for every i}. 

Note that Sq is cr-invariant. We will construct a continuous "almost one-to-one" map /i : So — > /C such 
that ho a = f o h. 

Claim. For each sequence a := (ai)igp} in So, there exists a unique point p G fC such that f^{p) G 
for every i >Q. 
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In order to prove the existence of p, one just needs to notice that the image under / of each of the 
intervals /i, I2, 13 is the union of the two other intervals. This implies that the intersection Di^o / ' i^ii ) 
is non-empty for every N, and therefore, that the intersection HieN f^^i^cn) is non-empty. The existence 
of p follows. In order to prove the uniqueness of p, observe that: for every e > 0, there exists v{s) > 1 
such that > ^{e) for every x G /C such that dist(x,Tj) > e for i = 1,2,3. Hence, if 

p p' were two points such that P{p) € /a, and P{p') £ for every i > 0, then one would have 
dist(/*(p), p{p')) —5- 00 ; this is absurd since the lengths of Ii, I2 and Is are finite. Hence there is at 
most one point p in HieN /~*(-^ai)- This completes the proof of the claim. 

Now, we consider the map /t : Sq ^ /C which maps a sequence a := (aj)igN to the unique point 
p G /C such that P{p) G /a^ for every i > 0. This map h obviously satisfies h o a = f o h. It is 
continuous (this is an immediate consequence of the continuity of /) and onto (because the image under 
/ of one of the intervals Ii, 12,13 is contained in the union of the two others). It is not one-to-one. 
For example, the Taub point T3 has two pre-images under h : the sequences (1,2,1,2,1,2,...) and 
(2,1,2,1,2,1,...). More generally, every point x £ IC such that (p) is a Taub point for some integer 
«o > (and thus /* (p) = /*° {p) for every i > io) has two pre-images under h, and these two pre-images 
are preperiodic for a. This is the only lack of injectivity of /i : if p G /C has not a single pre-image under 
h, then there exists io > such that /*" (p) is a Taub point (this follows from the fact that the intersection 
between two of the three intervals hjh, h is reduced to a Taub point). 

It follows that the image under h of a closed cr-invariant aperiodic subset of Sq is a closed forward- 
invariant aperiodic subset of /C. So we are left to prove that the union of all the closed a-invariant and 
aperiodic subsets of Sq is in dense in Eq. 

A element a = (ai)i>o of Sq is said to be square-free if it does not contain the same word re- 
peated twice : for every iq > and every £ > 0, the word ajg . . . ajg+^-i is different from the word 
aj,i+£ . . . ai(,+2^-i- It is well-known that there exist square-free elements in Eq (such an element may be 
easily deduced from the well-known Prouhet-Thue-Morse sequence, see for example [1, corollary 1]). 
Now let a = {ai)i>o G Eq be square-free, then the cr-orbit of a does not accumulates on any periodic a- 
orbit, hence the closure of the cr-orbit of a is a closed forward-invariant aperiodic subset of Eq. Moreover 
the same is true, if one replaces a by a sequence a' G Eq which has the same tail as a {i.e. there exists zq 
such that a'j = for i > io). The set of all sequences a' which has the same tail has a is obviously dense 
in Eq. Hence the union of all closed cr-invariant aperiodic subsets of Eq is dense in Eq. As explained 
above, the proposition follows. □ 
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Abstract. We study the asymptotic behavior of vacuum Bianchi type A spacetimes. It has been conjec- 
tured that this behavior is driven by a certain circle map, called the Kasner map. As a step towards this 
conjecture, we prove that some orbits of the Kasner map do indeed attract some solutions of the system of 
ODEs which describes the behavior of vacuum Bianchi type A spacetimes. The orbits of the Kasner map 
for which we can prove such a result are those which are not periodic and which do not accumulate on 
any periodic orbit. This shows the existence of Bianchi spacetimes with a aperiodic oscillatory asymptotic 
behavior. 



1. Introduction 

A Bianchi spacetime is a cosmological spacetime which is spatially homogeneous. More precisely, a 
Bianchi spacetime is a spacetime (M, g) such that M is diffeomorphic to the product G x / of a three- 
dimensional Lie group G and an interval / of R, and g = hf — dt^ where ht is a left-invariant riemannian 
metric on G x {t} ~ G for every tel. Hence a Bianchi spacetime can be seen as a one-parameter 
family {ht)tei of left-invariant riemannian metrics on a given three-dimensional Lie group G. A vacuum 
Bianchi spacetime is a Bianchi spacetime satisfying the vacuum Einstein equation Ric^ = 0. A Bianchi 
spacetime is of type A if the corresponding three dimensional Lie group G is unimodular (otherwise it is 
of type B). 

The space of left-invariant riemannian metrics on a given Lie group is finite-dimensional. Therefore, 
in the restricted context of Bianchi spacetimes, the Einstein equation is simply an ODE on a finite dimen- 
sional phase space. A classical and convenient way to write this ODE is to use the so-called Wainwright- 
Hsu variables (Si, E2, S3, A^i, A^2, ^^3). The corresponding phase space is the four-dimensional mani- 
fold 

H = {(Si, S2, S3, ATi, iV2, 7V3) G I Si + S2 + S3 = , n = 0} 

where 

n = 6 - (S? + si + si) + ^(ATf + Ar| + ATf) - + NiN^ + N2N^). 
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and the vacuum Einstein equation reads 



(S) 



1 


= (2-a)Si--Ri 




= (2-g)S2-i?2 




= (2-g)S3-i?3 


N[ 


= -{q + 2^i)Ni 




= -{q + 2^2)N2 




= -{q + 2E3)N3 



where 



q = 




T.I + E^) 








Ri = 


I 


- iV| - N'i 


+ 2N2NS 




- N1N2) 


R2 = 




- iV| - Nf 


+ 2A^3A^i 


- A^2iVi 


- N2N3) 


R3 = 




- Nf - Ni 


+ 2N1N2 


-N3N2 


- N3N1) 



(see, for example, [3] for the construction of the Wainwright-Hsu variables and the derivation of the 
above equations). In other words, vacuum Bianchi type A spacetimes can be seen as the parametrized 
curves in B that are solutions of the above system of ODEs. We will call Wainwright-Hsu vector field, 
and denote by X, the vector field on B corresponding to this system of ODE. We will denote by X* 
the time t map of the flow of X. We will denote B~^ the subset of B made of the points for which 
A^i, N2 and A''3 are non- negative. To measure the distances on B, we will use the riemannian metric 
h = 6,1:1 + dY^l + dS| + dNl + dN^ + diV| on B. 

Remark 1.1. We have chosen the "anti-physical" time orientation: Bianchi spacetimes will be future 
incomplete, but not necessarily past incomplete. With this "anti-physical" convention, the so-called 
mixmaster attractor (see below) is indeed an attractor (and not a repellor). 

The phase space B admits a natural stratification, according to the signs of the variable A'^i, N2, N3, 
this stratification being invariant under the flow of X. The different strata correspond to the different 
types of three-dimensional Lie groups. The behavior of the orbits of X contained in the one, two and 
three-dimensional strata (so-called, type I, II, VIq and VIIq orbits) is fully-understood. On the contrary, 
there is almost no rigorous result on the behavior of the orbits in the four dimensional biggest strata 
(co-called type VIII and IX orbits). 

The smallest stratum is called the Kasner circle, denoted by Bi or /C, and corresponds to the case 
where the Lie group G is M^. It is made of the points of B for which Ni = N2 = N^ = 0: 



/C = {(Si,S2,S3,7Vi,iV2,7V3) gH; TVi = 7V2 = 7V3 = 0} 

= {(Si,E2,S3,iVi,iV2,Ar3)GM6; iVi = iV2 = iV3 = Si + E2 + S3 = q = 2} 
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This is a euclidean circle in M^. The orbits of X contained in iSj = /C are called type I orbits. The vector 
field X vanishes on /C. Hence, every type I orbit is a fixed point. There are three points on the Kasner 
circle, called the special points or the Taub points which play a very important role in every attempt 
to understand the behavior of the solution of the Wainwright-Hsu equations; these are the points for 
which (El, S2, S3) is equal to (2, —1, —1) or (—1, 2, —1) or (—1, —1, 2); there are usually denoted by 
Ti, T3. For every point p G /C, the derivative DX{p) has four distinct eigendirections: the direction 
tangent to the Kasner circle which corresponds to a zero eigenvalue since X vanishes on /C, and the three 
directions and which are respectively associated to the eigenvalues —(2 + 2Ei), —(2 + 

2S2) and — (2 + 2S3). If p is one of the three special points Ti,T2, T3, then two of the three eigenvalues 
— (2 + 2S1), — (2 + 2S2), — (2 + 2S3); hence, the derivative of DX{p) has a triple zero eigenvalue. If pis 
not one of the three special points Ti, r2, T3, then the eigenvalues -(2 + 2Si), -(2 + 2E2), -(2 + 2E3) 
are pairwise distinct, two of them are negative, and the other one is positive; in other words, the derivative 
DX{p) has two distinct negative eigenvalue, a multiplicity one eigenvalue, and a positive eigenvalue. 

The next stratum is denoted by Bjj and correspond to the case where G is the Heisenberg group. It is 
made of the points of B for which exactly two of the ATj's are equal to zero: 

^iU% = {(Ei,E2,E3,A^i,iV2,7V3) A^i =iV2 = 0oriVi = Ar3 = 0orAr2 = iV3 = 0} 

One can easily check that Bi U 6jj is the union of three ellipsoids which intersect along the Kasner circle 
/C = The vector field X commutes with the transformation {Ni,N2, N3) i-> (— A^i, —N2, —N^), 
so we can work in Hjj fl to simplify the discussion. The intersection Hjj fl is the union of three 
disjoint hemi-ellipsoids bounded by the Kasner circle. The orbits in B2 are called type II orbits. Type 
II orbits can be calculated explicitly (see for example [3]). It appears that every type II is a heteroclinic 
orbit connecting a point of /C \ {Ti, T2, T3} to another point of /C \ {Ti, T2, r3}. Moreover, for point 
p G /C \ {Ti, T2, T3}, there is exactly one type II orbit in 6+ which "starts" at p (this orbit is tangent to 
the eigendirection associated to the positive eigenvalue of DX{p)) and two type II orbits in 6+ which 
"starts" at p (these orbits are tangent to the eigendirections associated to the negative eigenvalues of 
DX{p)). The Kasner map / : /C ^ /C is defined as follows. The three special points Ti, T2, T3 are fixed 
under /. For every point p G /C \ {Ti,T2,T^}, one considers the unique type II orbit which "starts" at 
p; this heteroclinic orbit connects p to another point q G /C \ {Ti, T2, T3}, which is by definition the 
point f{p). We will denote by O^ jf^p^ the unique type II heteroclinic orbit connecting p to f{p) in B^. 
The Kasner map / admit a nice geometric description. Consider the equilateral triangle which is tangent 
to K, at the three special points Ti,T2,T3. Denote by Mi,M2,M3 the vertices of this triangle. For 
p G /C \ {Ti, T2, T3}, let Mi be the vertex of /C which is the closest to p. The line (Mjp) and intersects 
the circle /C at two points: the point p and the point f{p). See figure 1. Using this geometric description, 
it is easy to see / is . Let us consider the derivative of / with the respect to the arc length on K, defined 
by the riemannian metric h. Then the norm of the derivative of the Kasner map at the three special points 
is equal to 1. But, for every p G /C \ {Ti, T2, r3} the modulus of the derivative / at p is strictly bigger 
than 1. The Kasner map / is non-uniformly expanding : for every compact set C in /C \ {Ti, T2, r3} 
there exists a constant f c > 1 such that || \Df{p) || |/j > z/c for every p in C. 

The strata B\ji^^ and -Byilo three-dimensional; Byy^^ is made of the points of B for which exactly 
one of the Ni& is equal to zero, the two others being of different signs; ^vilo made of the points 
of B for which exactly one of the A'j's is equal to zero, the two others being of the same sign. They 
correspond respectively to the case where G = 0(1, 1) ix and G = 0(2) ix M^. Observe that Hyig 
does not intersect B'^. 
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Figure 1 . Construction of an orbit of the Kasner map 

Finally the strata Bym and Bix are four-dimensional (these are open subsets of B). They correspond 
to the points of B for which none of the Ni 's is equal to zero. In ^yill '^^ hawe the same 

sign, whereas in Bix the three Ni's have the same sign. These strata correspond respectively to the cases 
where G = SL(2, M) and G = S0(3, M). Observe that ^yill does not intersect B+ 

Misner has conjectured that the dynamics of type IX orbits "is driven" by the Kasner map ([5]). The 
idea is that each "heteroclinic chain" (7UOq j(q) U/(g)UOj(g) j2(g) U. . . should attract some Bianchi IX, 
and that "generic" Bianchi IX orbits should correspond to "generic" orbits of the Kasner map (see [3] for 
a detailed discussion various possible precise meanings of this conjecture) (recall that by 
is the type II heteroclinic orbit connecting /"(g) to /"^^(g) in B^. 

In 2001, Ringstrom proved that A := JC U Bn is a global attractor : for "almost every" point q in 
Bix, "^he distance from the orbit of q to A tends to as the time goes to oo ([9], see also [4]). Although 
remarkable, Ringstrom's theorem does not tell anything on the relation between the dynamics of Bianchi 
orbits and the dynamics the Kasner map. The purpose of the present paper is to prove that "there are 
many point q £ IC such that the 'heteroclinic chain q U Oq, f{q) U f{q) U C/(g),/2(^) U . . . does attract 
some Bianchi IX orbits". 

Definition 1.2. For g G /C, we call stable manifold of q, and we denote by W^{q), the set of all points 
r € ;B for which there exists an increasing sequence of real numbers (tn)n>o such that 

dist(X*"(r),r(g)) ^ 

n— >oo 

and such that the Hausdorff distance between the piece of orbit {X*(r) ; tn <t < tn+i] and the type II 
heteroclinic orbit C'jn(g) jn+ij^g) tends to as n — oo. 

A subset C of /C is said to be f -invariant if /(C) C C. It is said to be aperiodic if it does not contain 
any periodic orbit for /. Our main result can be stated as follows: 

Tlieorem 1.3. If q G K, is contained in a closed f -invariant aperiodic subset of IC, then the stable 
manifold W^{q) contains a embedded three-dimensional disc Dj^{q) embedded in Bjx- This disc 
can be chosen to depend continuously on q when q ranges in a closed /-invariant aperiodic subset ofJC. 
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Remark 1.4. The analog of theorem 1.3 for Bym instead of Bix also holds. We will focus on Bix to 
simplify the notations. 

Observe that the hypothesis of theorem 1.3 is satisfied by a dense set of points in /C : 

Proposition 1.5. The union of all the closed f -invariant aperiodic subsets ofK, is in dense in /C. 

Remark 1.6. Let E be the union of all the closed, forward invariant, aperiodic subsets of K. According to 
the above proposition, E is dense in /C. Observe nevertheless that is a "small" subset of K, both from 
the topological viewpoint (it is a meager set) and from the measurable viewpoint (it has zero Lebesgue 
measure). Also observe that theorem 1.3 does not tell that W'^{q) depends continuously on q when q 
ranges in E. So, we do not know whether the union of all the stable manifolds W^{q), where q ranges 
E, is dense in B (or in an open subset of B) or not. 

The origin of the hypothesis of theorem 1.3 is purely technical: if p belongs to a /-invariant aperiodic 
compact subset of /C, we can find a coordinate system in which the vector field X is "linear"; this makes 
the estimates on the flow of X near such a point much easier. We do not know if such linearizing 
coordinate systems exist near a point of /C which is periodic or preperiodic under the Kasner map. 
Nevertheless, it should be noticed that the estimates on the flow of X that we need to construct stable 
manifolds are much weaker than the existence of a linearizing coordinate system. So, we do think that 
theorem 1.3 can be extended to any point g G /C such that the orbit of q under the Kasner map does 
not accumulate one any of the three special points Ti, T2, T3. Understanding the behaviour of the orbits 
which pass arbitrary close to the three special points seems to be a much harder problem. 

As we were finishing to write the present paper, M. Georgi, J. Harterich, S. Liebscher and K. Webster 
put on arXiv a preprint in which they prove that the period three orbits of / admit non-trivial stable 
manifolds ([2]). At the end of this preprint, they claim that their techniques can be used to extend their 
result to any orbit of / which does not accumulate one any of the three special points Ti, T2, Ts. Such a 
extension would imply our result. 

Let us try to give a rough idea of our proof of theorem 1.3. Consider a point p on the Kasner circle 

K, which is not preperiodic under the Kasner map /. Using the link between the behavior of the orbit of 
p under /, the continued fraction development of the Kasner parameter u{p), and the expression of the 
eigenvalues of DX{p) in terms of the Kasner parameter u{p), we show that the eigenvalues of DX{p) 
are independent over Q. Then, using a generalization of Sternberg's theorem due to Takens, we obtain 
a coordinate system on a neighborhood VpOfp such that X is "essentially linear" in this coordinate 
system. Since X is "essentially linear", we can compute explicitly the flow of X. This leads to the 
following important observation: when an orbit of X passes in the neighborhood Vp, the distance of this 
orbit to the mixmaster attractor decreases super-linearly, while the drift in the direction tangent to the 
attractor is as small as one wishes. Now consider a closed /-invariant aperiodic subset C of /C. One 
can find a finite number of points pi , . . . , p„ in C such that the neighborhoods , . . . , Vp„ cover C. 
When a type IX orbit goes through one of the neighborhoods l^Ji > • • • ) ^„ the distance of this orbit to the 
mixmaster attractor decreases super-linearly. Of course, this orbit will also do some excursions outside 
U • • • U Vp^, but the duration of these excursion is uniformly bounded (at least when the orbit is 
close to a type n heteroclinic orbit connecting a point of C to another point of C). So the effect of 
these excursions outside , . . . , Vp^ can be neglected when compared to the super-linear contraction 
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that occurs inside , . . . , Vp^ . All in all, this will allow us to show that the drift of the orbits in the 
direction tangential to the mixmaster attractor A can be neglected when compared to the contraction of 
the distance to this attractor. This will lead the existence of stable manifolds. 

Sections 2.. .7 are devoted to the proof of theorem 1.3 . Proposition 1.5 will be proved in section 8. 

Acknowledgements. I would like to thank Lars Andersson for some stimulating discussion, and for his 
encouragements to write the present paper. 

2. Takens linearization theorem 

Let y be a C°° vector field on a manifold M, and p be a point of M. such that Y{p) = 0. The linear 
space TpM admits a unique decomposition 

where Ep, Ep, Ep are -invariant linear subspaces, the real parts of the eigenvalues of the restric- 

tion of DY{p) to Ep (resp. E^, Ep) are negative (resp. equal to 0, positive). We shall denote by s, c and 
u the dimensions of the linear subspaces Ep, Ep and Ep. 

Definition 2.1. The vector field Y satisfies Sternberg's condition atp if the eigenvalues of DY {p)\E^^E^ 
(counted with multiphcities) are independent over Q. 

F. Takens as proved the following generalization of the classical Sternberg's linearization theorem : 

Theorem 2.2 (see [10], page 144). Assume that Y satisfies Sternberg's condition at p. Then, for every 
r > 0, one can find a neighborhood U and a C coordinate system (xi, . . . , iCs, yi, . . . , tjci -^i, • • • , Zy^ 
on U centered at p, such that, in this coordinate system, Y reads: 

^ d ^ d ^ d 

(1) Y=Y1 aij{yi, yc)xj^ + (t>i{yi, • • ■ > yc) o- + ^ h,j{yi, yc)zjg- 

t,J=l 1=1 lj=l 

where the eigenvalues of the matrix (aij(0, . . . , 0))i<jj<s have negative real parts, the eigenval- 
ues of the matrix {d4>i/dyj{0, . . . ,0))i<ij<c o.re on the imaginary axis, and the eigenvalues of 
{bij{0, . . . , 0))i<jj<„ have positive real parts. 

Remarks 2.3. In general, the size of the neighborhood U does depend on the integer r in general, and it 
is not possible to find any C°° local coordinate system centered at p such that (1) holds. 

The name "Takens linearization theorem" is shghtly incorrect: the vector field Y is not linear in 
the {xi, ... ,Xs,yi, ■■■ ,yc, zi, Zu) coordinate system. Nevertheless, observe that the submanifold 
defined by (j/i = • • • = j/c = 0) is invariant under the flow of Y, and the restriction of Y to this subman- 
ifold is linear. For (yj', . . . , y[?) 7^ (0, . . . , 0), the submanifold defined by (yi, . . . , yc) = (y?, . . . , y^) 
is not invariant under the flow of Y in general, but the projection of Y on the tangent space of this 
submanifold is linear. 

The vectors ^(p), ■ ■ ■ , g|"(p) span the subspace Ep, the vectors -^(p), • • • > ^(p) span the sub- 
space Ep and the vectors -^{p), . . . , ^(p) span the subspace Ep-. this follows from equality (1) 
together with the signs of the real parts of the eigenvalues of the matrices (aij(0, . . . , 0))i<jj<s, 
(90j/%(O, . . . , 0))i<jj<c and (6ij(0, . . . , 0))i<i,j<u. 



APERIODIC OSCE^LATORY ASYMPTOTIC BEHAVIOR FOR SOME BIANCHI SPACETIMES 



7 



3. Linearization of the vector field X near a point of the Kasner circle which is 

NOT PREPERIODIC 

Recall that X denotes the vector field on B whose orbits are the solution of the system of ODEs (S). 
We will prove that the vector field X satisfies Sternberg's condition at every point of the Kasner circle 
which is not preperiodic under the Kasner map /. For this purpose, we need to introduce the so-called 
Kasner parameter. \f q = S2(g), S3(g), 0, 0, 0) is a point of the Kasner circle, the Kasner 

parameter u{q) G [1, oo] of q is defined by 

O) 2-EU.;)s.WS3W = 27 ,"y<i + «fa)); 

(1 + u{q) + u{qY)'^ 

The map q >->■ is not one-to-one: indeed the left hand side of the above equality is invariant under 
permutation of the Tii{q)^ ^2{q)j ^3(5). Nevertheless, ^2{q)j S3 ((7) are characterized up to per- 

mutation by u{q). Hence, the eigenvalues -(2 -|- 2Ei(g)), -(2 + 2E2(g)), -(2 -1-2X3(9)) of DX(g) are 
characterized up to permutation by the Kasner parameter u{q). More precisely: 

(3) { - (2 + 2Si(g)), -(2 + 2S2(g)), -(2 + 2X3(5))} 

u{(i) l + u{q) u{q){l + n{q)) 



l+u{q)+v?{q)' l+u{q)+v?{q)' l+u{q)+v?{q) 

(this can be deduced from (2), together with the equalities ^i{q) + ^lifl) + ^3(5) = and '^\{q) + 
Ti^iq) + ^Kg) = 6). Note that u{q) = 00 if and only if q is one of the three special points Ti, T2, T3. 
The main advantage of the Kasner parameter is the fact that the Kasner map / admits a nice expression 
in terms of this parameter. For every g G /C, one has u{f{q)) = / {u{q)) where / : [1, +00] — > [1, 00] is 
defined by 



(4) Ru) = { 



00 if -u = 1 or 00 
u-l if u{q)> 2 
^ «^ if !<«(?)< 2 



From now on until the end of the section, we consider a point p G /C \ {Ti, r2, T3}. As we already 
pointed out in the introduction, the derivative DX{p) : TpB TpB has a four distinct (multiplicity 
one) real eigenvalues: two negative eigenvalues (corresponding to the directions of the two type II in 
B~^ which "arrive" at p), a multiplicity one zero eigenvalue corresponding to the direction of the Kasner 
circle, and one positive eigenvalue corresponding to the unique type n orbit in B~^ which "starts" at p. 
We will denote the non-zero eigenvalues of p by Xsip), A'^(p), Au(p) with Xs{p) < Xg{p) < < Xu{p)- 
Recall that (A,(p), A',(p), A„(p)) = ((2 + 2Si(g)), -(2 + 2X2(5)), -(2 + 2X3(5))) up to permutation). 
From (2), we easily deduce: 

c^^ f \ ( \ \' ( \ \ ( \ \ - ( + 1 + u{p) u{p) \ 

PJ yAsm . \KP) , ^uKP) ) 1^ i + ^(^)^^2(p)' l + u{p)+U-'{p)' 1 + U{p)+V?{p)) 

The following observation will play a key role in the proof of our main theorem 1.3: 
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Proposition 3.1. The contraction rate of DX{p)\E^ dominates the dilatation rate of DX{p)\E^, that is 



Up) 



> 1 and 



Kip) 



K{p) 



> 1. 



Proof. This is an immediate consequence of (5) and of the fact that > 1 and "^"""^"^ = 1 + u > 1 
for every u > 1. □ 

The statement below gives some necessary and sufficient condition for the vector field X to satisfy 
Sternberg's condition at p, in terms of the behavior of the orbit of p under the Kasner map. The hypothesis 
of our main theorem 1 .3 comes directly from this condition. 

Proposition 3.2. The three following conditions are equivalent : 

(1) the vector field X satisfies Sternberg's condition at p ; 

(2) the Kasner parameter u{p) is neither a rational, nor a quadratic irrational number ; 

(3) p is not a preperiodic point for the Kasner map f. 



Proof. Let us prove I 4^ 2. The vector field X satisfies Sternberg's condition at p if and only if the 
eigenvalues Xs{p), Kip): Xu{p) are independent over Q. Using formulae (5), one sees that this is equiv- 
alent to the fact that the numbers u{p), 1 + u{p) and u{p){l + u{p)) are independent over Q. Clearly, 
this is equivalent to the fact u{p) is neither a rational number nor a quadratic irrational number. 

Now, let us prove 3 4^ 2. For this purpose, we will use the action of Kasner map / on the Kasner map. 
Recall that, for every q on the Kasner circle, one has u{f{q)) = f {u{q)) where / : [1, +oo] [1, oo] is 
given by (4). Every orbit of / which does not end on oo, passes an infinite number of times in the interval 
[1, 2]. Let k : (1, 2] (1,2] be the return time function of / in the interval (1, 2], and F : (1, 2] (1, 2] 
be the first return map of / in the interval (1,2], that is 



k{u) = inf{n > such that /"(ti) G (1, 2]} and F{u) = f^^'^\u). 



Then, 



k{u) 



u-l 



- 1 and F{u) 



u-l 



+ 1, 



where [xj is the integer part of x and {x} := x — [xj is the fractional part of x. Observe that a point 
p G /C is preperiodic for the Kasner map / if and only if u{p) is a preperiodic point for F. Observe that 
F is just the Gauss map u i-^ conjugated by the translation u u + 1, and that k{u — 1) + 1 
is the first term of the continuous development fraction of u. On the one hand, it is well-known (and 
very easy to prove) the preperiodic points of the Gauss map ?i { } are exactly the real numbers 
u G (0, 1] such that the sequence of integers which appear in the continuous fraction development of u is 
finite or preperiodic. It follows that a point p G /C is preperiodic for the Kasner map / if and only if the 
continuous fraction development of u{p) is finite or preperiodic. On the other hand, it is well-known that 
the continuous fraction development of u G M is finite if and only if u is a rational number, and that the 
continuous fraction development of u is preperiodic if and only if it is a quadratic irrational number. This 
shows that p G /C is a preperiodic for the Kasner map / if and only if u{p) is a rational or a quadratic 
irrational number. □ 
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According to proposition 3.2, if p is not a preperiodic point for the Kasner map, the hypotheses of 
Takens linearization theorem are satisfied by the vector field X at p. This theorem provides us with local 
coordinates near p in which the vector field X has a very simple expression : 

Proposition 3.3. Assume that the point p is not eventually periodic under the Kasner map. Then there 
exists a neighborhood U ofp in B and a coordinate system {x, x' , y, z) on U centered p, centered at 
p, and such that, in this coordinate system, X reads : 

d d 3 

(6) X{x, x', y, z) = Xs{y)x— + X'g{y)x'— + Xu{y)z— 

where Xs{y) < Kiv) < ^ '^"'^ ^uiv) > ^far every y. 

Remarks 3.4. Let p, U, {x, x', y, z), Xs{y) , X'g{y) , Xu{y) be as in the above proposition. 

(i) The vector field X vanishes on the one-dimensional submanifold {x = x' = z = 0} and 
nowhere else. It follows that this submanifold is the intersection of the Kasner circle with U : 

jcnu = {x = x' = z = o}. 

(ii) For every y, the real numbers Xs{y) , X'g{y) , X^iy) are the three non-zero eigenvalues of the de- 
rivative of X at the point of coordinates (0, 0, y, 0) € K,. Recall that this derivative also has one 
zero eigenvalue (the corresponding eigendirection being tangent to the Kasner circle). 

(iii) For every c, the three-dimensional sub-manifold {y = c} is invariant under the flow of X and 
the restriction of X to this submanifold is linear. 

(iv) The three sets {x = 0}, {x' = 0} and {z = 0} are three-dimensional submanifolds, which 
are invariant under the flow of X, and contain JCflU. It follows that, up to permutation, these 
three submanifolds correspond to the submanifolds {Ni = 0}, {N2 = 0} and {N3 = 0}. In 
particular, 

(HvHo U % U /C) n [/ = {a; = 0} U {x' = 0} U {z = 0}. 

(v) As a consequence, the three two-dimensional submanifolds {x = x' = 0}, {x = z = 0} 
and {x' = z = 0} coincide, up to permutation, with the three two-dimensional submanifolds 

{Ni = N2 = 0}, {iVi = iVs = 0} and {iV2 = = 0}. In particular, 

(% n /c) n ;7 = {z = x = 0} u {z = z = 0} u {z' = z = 0}. 

(vi) The vector field Xs{y)xd / dx + X'g{y)x'd/dx' + Xu{y)zd/dz is not changed if we multiply the 
coordinates x',y,z by some arbitrary real numbers. Together with item (iv), this shows that, 
up to multiplying x and/or x' and/or z by —1, we may assume that 

B+nU = {x>Q,x' >0, z>0}. 

(vii) On the one hand, since K. H U = {x = x' = z = 0}, the metric induced on /C fl C/ by the 
riemannian metric dx^ + dx^ + dy"^ + dz'^ is simply dy"^. On the other hand, the vector field 
Xs{y)xd / dx + X'g{y)x'd/dx' + Xu{y)zd/dz is not changed if we replace the coordinate y by 
ip{y) where : M — > M is any C°° diffeomorphism such that (p{0) = 0. Therefore, up to 
replacing y by ip{y) for some diffeomorphism (p, we may assume that the metrics induced on 
JCr\Uhydx'^ + dx'"^ + dy"^ + dz^ and h = dT,l + + c?S| + dNf + dN^ + dN^ coincide. 

Proof of proposition 3.3. The derivative DX(jp) has two negative, one zero, and one positive eigenvalue. 
According to proposition 3.2, the vector field X satisfies the Sternberg condition at^?. Therefore, a crude 
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application of Takens' theorem 2.2 implies that there exists a local coordinate system {xi,X2,y, z) 
on a neighborhood U oip 'va B, centered at p, such that : 

2 2 Q Q Q 

(V) x{xi,x2,y,z) = ^Yl «r: + '^(2/) o- + Ky)z-Q- 

for some real valued functions ai^i, ai 2, 02,1, «2,2! ^' defined on a neighborhood of in M. Moreover, 
the eigenvalues of the matrix (aij (O)) are negative, and 6(0) is positive. Replacing U by & smaller 
neighborhood of p if necessary, we can assume that the three Taub points Ti,T2, are not in U. 

Now, (7) implies that the curve {xi = X2 = = 0} is the only curve in U containing the point p, 
invariant under the flow of X, and such that DX{p) vanishes on the tangent space at p of this curve. 
Hence, the curve {xi = X2 = z = 0} has to be the intersection of the Kasner circle /C with U. Since X 
vanishes on /C, it follows that = 0. 

For j/° small enough, let qyO be the point of coordinates {xi,X2,y, z) = (0, 0, ?/°, 0). This is a point of 
the Kasner circle /C, which is not a Taub point. Hence, the point QyO has four distinct eigenvalues : two 
distinctive negative eigenvalues \s{y^) < Kiv^) < 0' 01^^ ^^I'o eigenvalue whose associated eigendi- 
rection is tangent to the Kasner circle, and one positive eigenvalue Xu{y^) > 0. The set {y = y^) 
is a three-dimensional manifold, transversal to the Kasner circle. Looking at (7), we see that this 3- 
dimensional submanifold is invariant under the flow of X, and that the restriction of X to this invariant 
manifold is linear.This shows that 6(c) = Xu{y^), and that Xs{y^) and A'^(y°) are the eigenvalues of the 
matrix (aij(y°))ij=i,2- Since \s{y^) and X's{y^) are distinct, there exists a linear change of coordinates 
{xi,X2, z) — > (x, x' , z) on the submanifold {y = y°), so that 

X{x,x',y',z) = X,{y^)x-^ + X'^{y^)x'-^+X^{y^)z^^. 

Since eigenvalues and eigendirections of the point QyO depend in a smooth way of c, one may per- 
form the above change of coordinates simultaneously for every y^, and get a system of C°° coordinates 
(x,x',y,z) defined on U, such that 

c) d d 

X{x, x', y, z) = Xs{y)x— + X'g{y)x'— + A„(y)z— . 

The proposition is proven. □ 

From now on, for every point p on the Kasner circle /C which is not preperiodic under the Kasner map, 
we fix a local coordinate system {xp, x'p, yp, Zp) on a neighborhood Up of p in B, centered at p, and 
such that, in this coordinate system, X reads : 

d d d 

(8) X{xp, x'p, yp, Zp) = Xs{yp)xp— + K{yp)Xp— + Xu{yp)zp— 

with Xs{yp) < X'g{yp) < and Xu{yp) > for every yp. According to items (i), (ii), (iv) and (v) of 
remarks 3.4, we automatically have 

(9) JCnUp = {x = x' = z = 0} 

(10) (%U/C)ni7p = {x = x' = 0}U{x = z = 0}U{x' = z = 0}, 

(11) (Hvii, u u /C) n [/p = {x = 0}u{x' = 0}u{z = 0}, 
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and XsiUp), KiVp): ^uiUp) are the three non-zero eigenvalues of the point (0, 0, yp, 0) G /C. Moreover, 
according to items (vi) of remarks 3.4, we may — and we will — assume that 

(12) B'^ nUp = {xp > , x'p > , Zp > 0}. 
We will consider the riemannian metric gp on Up defined by 

(13) gp := dxp + dxp + dyp + dzp. 

Recall that h is the riemannian metric on B defined by h = dSf + dSl + dSg + dN^ + d7V| + dN^. 
According to item (vii) of remarks 3.4, we may — and we will — assume that 

(14) the riemannian metrics induced by gp and h on ICn Up coincide. 
For a < P and 7 > 0, we consider the subset Vp{a, /3, 7) of Up defined as follows: 

(15) yp(a,/3,7) := Up n{0< Xp,Xp,Zp < 7 and a<yp< /S}. 

Observe that Vp is contained in Up n B^. If a, (3, 7 are small enough, then Vp{a, 7) is contained in the 
interior of Up. \i a < d < (3, then Vp{a, /3, 7) is a neighborhood of p in B~^. 

4. DULAC MAP NEAR A POINT OF /C WHICH IS NOT PREPERIODIC 

Let p be a point of the Kasner circle /C which is not preperiodic for the Kasner map. We consider the 
neighborhood Up of p and the local coordinate system {xp, Xp, yp, Zp) on Up defined at the end of the 
preceeding section. In particular, we assume that properties (8). ..(12) and (14) are satisfied. Using (8), 
one can calculate explicitly the time t map of the flow of the Wainwright-Hsu vector field X in the 

(xp, x'p, yp, Zp) coordinate system : 

(16) X (^Xp, Xp, yp, Zp) = (^XpC "^^p^ , Xp€. ^^^^^ , yp , ZpC "(^f^ ^ . 

Now, we consider some real numbers a < (3 and 7 > small enough so that the set 

(17) Vp := Vp{a, (3, 7) = {0 < Xp, x'p, Zp < j and a < yp < /?} 

is contained in the interior of Up. The purpose of this section is to analyze quite in details the behavior 
of the orbits of the flow of X inside Vp. In the [xp, x'p, yp, Zp) coordinate system, 1^ is a 4-dimensional 
box, i.e. the cartesian product of four closed intervals. The boundary of Vp is made of eight faces. Three 
of these eight faces will play an important role in the remainder of the paper: 

(18) Mp := n {xp = 7} M'p■.= Vpf^{x'p = -f} Np := Vp ^ {zp = -f} . 

Looking at (8), we notice that X is transversal to Mp, Mp and Np, and is tangent to the five other faces 
of Vp. Moreover, we notice that X is pointing inward Vp along Mp and Mp, it is pointing outward Vp 
along Np. It follows that: 

• an orbit of X can enter in Vp by crossing either the face Mp or by crossing the face Mp, 

• an orbit of X can only exit Vp by crossing the face Np. 

We need to analyze quite precisely the behavior of the orbits of X inside Vp. First, we need to 
understand which type n orbit intersect the hypersurfaces (with boundary and comers) Mp, Mp and Np. 
Recall that, for every point q G Hjj, the orbit of g is a heteroclinic orbit connecting two points of /C. 

Proposition 4.1. Let q be a point Bjj, and denote by O the orbit of q. Let a{q) = a{0) denote the 
unique a-limit point ofO, and oj{q) = co{0) denote the unique uj-limit point ofO. 
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(1) The orbit O intersects the hypersurface Np if and only if the point a{q) is in Vp. 

(2) The orbit ofO intersects the hypersurface Mp U Mp if and only if the point a;(g) is in Vp. 

Proof. We prove the first statement; the second one follows from similar arguments. Of course, we will 
work in the {xp, Xp, yp, Zp) coordinate system. According to (10), (15) and (18), 

/Cnl/p = {xp = Xp = 0, a < 2/p < = 0}, 
%nA'p = {xp = x'p = Q,a<yp<^, Zp = -f}. 

Suppose that the orbit O intersects the hypersurface Np at some point q = (0, 0, y, 7), with a < y < p. 
Then, according to (16), the past orbit of q is contained in Vp, and converges to the point a{q) = 
(0, 0, y, 0) G /C n Vp. In particular, the point a{q) is in Vp. Conversely, suppose that the point a{q) is in 
Vp. Then a{q) = (0, 0, y, 0) for some y G [a, /?]. Using again (16), we see that the only orbit of X in 
converging towards the point (0, 0, y, 0) as t — > —00 is the curve t ^ (0, 0, y, e^"^^')*). Hence, the orbit 
O intersects the hypersurface Np at the point q = (0, 0, y, 7). □ 

This proposition allows to define two maps 

ap-. B^nNp — > icnVp Up-. ^+ n {Mp u m^) — > icnVp 

q I — > a{q) q 1 — > u>{q) 

The map ap is one-to-one (there is only one type II orbit in B'^ which "starts" at a given point of /C), 
whereas the map ujp is two-to-one (there are two type II orbits in i3+ which "arrive" at a given point of 
/C). The restriction of Up to fi Mp (resp. to i3jj n Mp) is one-to-one. 



Proposition 4.2, The maps ap and LOp are local isometrics with respect to the metrics induced on 



Bjj n Np, Bjj n {Mp U Mp) and /C fl V^ by the riemannian metric gp = dXp + dXp + dyp + dZp 



Proof. The proof of proposition 4.1 shows that, in the {xp,Xp,yp, Zp) coordinate system, the map ap 
reads ap{0, 0, yp, 7) = (0, 0, yp, 0). Similarly, the map Up reads ojp{'y, 0, yp, 0) = (0, 0, yp, 0). □ 

Now, we have to define and study the Dulac map associated to Vp-. this map will describe the effect of 
going through Vp for the orbits of X. 

Let g be a point in Mp. If the Zp{q) > (which is the case if g G Byx), then (16) shows that the 
forward orbit of q will eventually exit Vp, by crossing the face Np. If Zp{q) = (which is the case if 
q G Bii), then (16) shows that the forward orbit of q will remain in Vp forever. Moreover, the orbit 
of q will converge towards the point u{q) G /C Pi V^. According to proposition 4.1, the type II orbit 
wiU eventually exit Vp, by crossing Np (since Q;(Oa;(g),/(u;(q))) = ^(9) ^ So, we may 

define a map 

$p : Mp — >Np 

as follows : 

• if Zp{q) > then ^p{q) is the first intersection point of the orbit of q with the hypersurfaceAfp ; 

• if Zp{q) = then ^p{q) is the first (and unique) intersection point of the type n heterocUnic orbit 
Ooj{q),f{oj{q)) with the hypersurface Np. 
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We call $p a Dulac map since it is the exact analog, in our situation, of the classical Dulac maps 
used to study planar vector fields. The above discussion together with formula (16) show that, in the 
(xp, x'p,yp,Zp) coordinate system, the map $p reads : 

(19) %{l,x'p,yp,zp) = 

(20) %{l,x'p,yp,0) = (0,0,2/p,l) 



^/ ^-A^feVA^y.) ^ ^ ifz„>0 



^p 



For every q G Mp, we decompose TqMp as a direct sum of two linear subspaces © where 



(21) 



(22) 



_d 
' dx: 



' dz. 



F^ := 



d 

'dyp 



iQ)- 



Similarly, for every r G Np, we write TrNp as a direct sum of two linear subspaces Gf. © where 



(23) 
(24) 



Gil 



d 
dxp 
d 

dyp 



(r) ©: 



(r). 



d 
'dx'p 



We are now ready to state the properties of the Dulac map $p which are the core of the proof of our 
main theorem 1.3: 

Proposition 4.3. The Dulac map <I>p : Mp — > Np is C^. Moreover, for every point q G Mp n {zp = 0}, 

the derivative D^p{q) : TqMp — )• T^^^^^-^Np of the map f&p at q satisfies: 

• D^p{q).v = for every vector V € F^; 

• D^p{q) maps F^ on G^^^^y and \\D^p{q).v\\gj^ = \\v\\gp for every vector v G F^. 

This proposition roughly says the following: when an orbit of X passes close to the point p G /C, 
the distance from this orbit to the mixmaster attractor A = IC U Bn is contracted super-linearly (in 
the {xp, x'p, yp, Zp) coordinates system, this distance in measured by the Xp, Xp, Zp coordinates), whereas 
there is no drift in the direction tangent to the attractor (in the (xp, Xp, yp, Zp) coordinates system, this 
drift is measured by the j/p coordinate). 



Proof of proposition 4.3. The fact that $p is (and even analytical) in restriction to Mp n {zp > 0} is 
an immediate consequence of formula (19). The fact that $p is on Mp n {zp = 0} follows from (19) 
and (20), together with the inequalities 

-M4>i and -#4>1 

(see proposition 3.1)^ The same ingredients actually show that, for every point q G Mp n {zp = 0}, 

d 



' dZn 



'dyp 



dyj 



^We do not know if $p is C^"*"^ for any e > 0, unless we have some better estimates on the ratios ^4!H and ^^i^ 



>'u(yp) 
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□ 



The coordinate Xp and Xp play similar roles. So, we will also consider the map 

: m; ^ Np 

defined as follows : 

• if Zp{q) > then ^p{q) is the first intersection point of the orbit of q with the hypersurfaceA/p ; 

• if Zp{q) = then ^'p{q) is the first (and unique) intersection point of the type II heterochnic orbit 

Oi^{q),f{oj{q)) with the hypersurface Np. 

In the {xp, Xp,yp, Zp) coordinate system, the map $p reads : 

(25) (x, , 1 , , Zp) = (,x,.z,-^=(^-)/^"(^'') , ,-K^y.)M,) , , i) if > 

(26) l,yp,0) = (0,0,yp,l) 

For every q G M'^, we will write TqM'p 

(27) F| ■-- 

(28) :-- 

Then, we can summarize the key properties of the map $p as follow : 

Proposition 4.4. The map $p : Mp — > Np is C^. Moreover, for every q G Mp n {zp = 0}, the derivative 
ofD^p{q) : TqMp T^i^^^q^Np of the map $p at q satisfies : 

• D^p(q).v = Ofor every vector v E ; 

• D^p{q) maps F^ on G^, and \\D^p{q).v\\gj^ = \\v\\gpfor every vector v in F^. 



5. Construction of a local section for a closed /-invariant aperiodic set of IC 

From now on until the end of section 7, we consider a closed, f -invariant, aperiodic subset C of the 
Kasner circle K. 

Observe that C is necessarly totally discontinuous. Indeed the points of /C that are preperiodic for the 
Kasner map / are dense in /C: this can be proved in several different ways; this follows for example the 
equivalence 2 3 of proposition 3.2. Let C be the union of C and of all the type II orbits connecting 
two points of C: 

qec 



as a direct sum of two linear subspaces F^ © F^ where 

d 

dyp 
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We want to prove that, for every point q E C, the stable manifold (q) contains a 3-dimensional disc 
(see definition 1.2 and theorem 1.3). To this end, we will consider a kind of "Poincare section" for C, 
and study the Poincare map associated with this section^. 

For every p G C, we consider the neighborhood Up of p, and the local coordinate system 
{xp, x'p, yp, Zp) on Up defined in section 3. Recall that, for every a, /3, 7, we denote 

Vpioi, P, 7) := Up n{0< Xp, Xp, yp<-faada<yp< P}. 

Recall that Vp{a, (3, 7) is a neighborhood of p in provided that a < < /3 and 7 > 0. For each 
p e C,we choose < < /3p and 7^ > small enough, so that 

i;0:=Fp(a°,/3°,7°) 

is contained in the interior of Up. Up to slightly modifying a^, /3p, 7p we can assume that the bound- 
ary of Vp is disjoint from C (that is, the points of coordinates {xp, Xp, yp, Zp ) = (0,0,q;,0) and 
(xp, Xp, yp, Zp) = (0, 0, /3, 0) are not in C): this is possible since C is totally discontinuous. 

Since C is compact, one can find a finite number of points pi, . . . ,pn & C such that the neighborhoods 
Vp^,. . . , Vp^ cover C. Now, we modify these neighborhoods in order to make them pairwise disjoint: 

• we set (api , f3p^ , 7pi ) := ("pi , Pp^ , 7pi )> and Vp^ := Vp, {ap, , pp, , 7pi ) = ; 

• then, we can find some constants a2,P2,j2 such that Vp,^ := Vp^{a2, ^2,12) is contained in 
Vp^ \ Vp^ , and such that C fl {Vp^ \ J is contained in the interior of Vp^ ; 

• then, we can find some constants a3,/33,73 such that Vp.^ := Ipjlaa, Asi Ts) is contained in 
Vp^ \ (Vp-, UVp^), and such that C n (1^° \ [Vp, U Vp^)) is contained in the interior of Vp^; 

• etc. 

At the end of this process, we get n pairwise disjoint domains Vp,, . . . , Vp^ , such that C is contained in 
the interior of U • • • U Vp^. For each i, there are some constants ai,(3i, 7, such that Vp^ = V^^'^^'"'\ 
and Vp^ is contained in the interior of Up^ . Hence, the result of section 4 apply to Vp^ . It may happen that, 
for some i, the point pi is not in Vp^ {i.e. that or /3j is non-positive), but we do not care. 

Now, we denote 

V:=Vp,U-- -UVp^ M:=Mp,U-. •UMp„ M' := M^^ U- ■ -UM^^ N := Np,U- ■ -UNp^. 

Then y is a neighborhood of C in B^. The hypersurfaces M, M' and N are transverse to X. An orbit 
of X can only enter in V by crossing M U M', and can only exit V by crossing A''. Moreover, according 
to proposition 4.1, we have the following important properties : 

Proposition 5.1. Every type II orbit whose u-limit point is in C intersects M U M'. Every type II orbit 
whose a-limit point is in C intersects N. 

We will see M U M' as a kind of "Poincare section" for C. Let us define the "Poincare map" Q 
associated to this section. First, we consider the "Dulac map" 

^■.MUM'^N 



The set C cannot admit a true Poincare section, since it contains some singularities of X (namely, the points of C). 
Nevertheless, we will consider a hypersurface with boundary N so that every type 11 orbit in C intersects N transversally. The 
hypersurface N will play the role of a Poincare section. 
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defined by ^\Mp^ = and ^\m^. = %i- For g 6 M U M': 

• if the forward of orbit of q exits V by crossing (which is the case if q € ^ix)^ ^^^^ ^il) 
definition the first intersection point of the orbit of q with the hypersurface A'^ ; 

• if the forward orbit of q remains in V forever (which is the case for every q e Bn), then is 
the first intersection point of the type II heterocUnic orbit C^j^^) with the hypersurface N. 

Now, we consider the "transition map" 

:N ^MUM' 

partially defined as follows: if the forward orbit of q re-enters in F by crossing the hypersurface MUM', 
then ^'(g) is the first intersection point of this forward orbit of q with M U M'; otherwise '^{q) is not 
defined. The "Poincare map" 

e-. MUM' ^ MUM' 
is by definition the product of the "Dulac map" $ and the "transition map" 

In the next section, we will study the dynamics Poincare map 0. For this purpose, we will use a 
riemannian metric on on B such that, for i = 1 , . . . , n, 

m = 9vi = dxl. + dx'l + dyl. + d^- 

6. Existence of stable manifolds for the Poincare map 9 : M u M' ->• M u M' 

The purpose of this section is to prove that, for every point g G C H (M U M'), the stable manifold of q 
for the "Poincare map" G contains a two-dimensional disc in M U M' . To this end, we will prove that 
C n (M U M') is a hyperbolic set for the map 0, and we will use a classical results on stable manifolds 
for hyperbolic sets. 

Definition 6.1. Let {M. , 5) be a riemannian manifold and (f) : M. ^ M. & map. A hyperbolic 
set for the map is a compact (/)-invariant subset C of such that, for every q ^ C, there is splitting 
TqM. = Fq ® Fq which depends continuously on q and such that, for some constant /i, G (0, 1) : 

(29) D(p{q).F^ C F^^^^ and \\D(j){x).v\\ < n\\v\\ for every qeCmdve F^ 

(30) D^{q).F^ = and \\{D^{q))-'^ .v\\ < iy\\v\\ for every qeC saidve F^. 

The dimension of the vector space F^ is called the index of C. The constant // is called a contraction 

rate of (p on C. 

Theorem 6.2. (see e.g. [6, page 167]) Let (j) : M ^ M be a map of a manifold A4, and C be a 
compact subset of M. which is a hyperbolic of index sfor the map (p. Then, for every e small enough, for 
every x E C, the set 

W^{x) ■.= {yeM\ dist{(j)^{y),(j)^{x)) < e for every n> 0} 

is a C^-disc of dimension s embedded in M, tangent to F^ at x, depending continuously on x (in the 
topology). Moreover, if fx is a contraction constant for (p on C, then there exists a constant k such that, 
for every e small enough, for every q £ C and every r G W^{q), 

distg{rir),riQ))<'^i^''- 
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We want to apply this theorem to the "Poincare map" 9 : M U M' M\J M'. So we need to prove 
that C n (M U M') is a hyperboUc set for 9. Recall that MUM' = Mp^ U • • • U Mp^ U M^^ U • • • U M^^ . 
So, for every q £ MU M', we have already defined a splitting TgM U M' = ® in section 4. It 
remains to prove that 9 satisfies (29) and (30) with respect to these splitting. For this purpose, we will 
use the decomposition of 9 as a product : 

9 = *o$. 

The behavior of the derivative of "Dulac map" $ was already studied in section 4 ; more precisely, we 
can rephrase propositions 4.3 and 4.4 as follows : 

Proposition 6.3. The map ^ : M U M' ^ N is C^. Moreover, for every q e C (1 {M U M'), the 
derivative ofD^{q) : Tq{M U M') T^^^-^N of the map ^ at q satisfies : 

• D^{q).v = Ofor every vector v G ; 

• D^(q) maps F^ on G^(gy cind \\^^{Q)-'^\\gp — ll'^llsp/'"' ^^^T vector v in F^. 

It remains to study the behavior of the derivative of the "transition map" ^ : iV — >• M U M'. We recall 
that ^{q) is well-defined only if the forward orbit of q intersects M U M'. So our first task is to show 
that * is well-defined at least on a neighborhood ofCCiNinN. 

Proposition 6.4. There exists a neighborhood VofCCiNin N, such that, for every q G V, the orbit 
of q intersects M U M' after some time t{q) which depends in a way on q. The map ^ is well- 
defined and on V. Moreover, there exists v > 1 such that, for every q G C Ci N, the derivative 
L>*(g) : TgN T^^q)M satisfies 

• D'^{q).Gq = F^^^^ and \\D'^{q).v\\g > v\\v\\g for every v G G^. 

Proof. Consider a point q e C H N. By proposition 5.1, the orbit of q intersects M U M' at some point 
r G C n (M U M'). Now, recall that : 

• N and M U M' are hypersurfaces with boundary that are transversal to the orbits of X ; 

• V was chosen so that C is contained in the interior of V. This implies that C does not intersect 
neither the boundary of the hypersurface A^, nor the boundary of hypersurface M U M'. Hence, 
q is in the interior of N, and r is in the interior of M U M'. 

This implies the existence of a neighborhood Vq of qin N such that, for every q' G V, the forward orbit 
of q' intersects M U M' after some time t{q') which depends in a way on q'. By definition of for 
every q' € Vq, we have ^(g') := X^^'^ Hq')- In particular, * is well-defined and on Vq. This proves 
the two first statements of the proposition. 

Since Bn is invariant under the flow of X, the map maps Bn n on Bn H (M U M'). Now, observe 
that, for every q € BuH N, the direction Gq is nothing but the tangent space of Bn (1 N at q, and the 
direction F^^^^ is nothing but the tangent space of Bji fl (M U M') at ^{q). This shows that d^{q) maps 
on F^(^) for every qeCHN. 

We are left tov > 1 such that > iy\\v\\g for every for every q e CCiN and every v G Gq. 

For this purpose, we will use the maps 

a:BiinN — >]CnV and u-.B^nN — >]CnV. 



18 



FRANCOIS BEGUIN 



We recall that a maps a point r £ B-^nN to the a-Umit point of the orbit of r, and that lo maps a point 
s G B]i n (M U M') to the w-Umit point of the orbit of s (see section 4). We also recall that a is a 
local isometry for the metrics induced by g on H N and ICOV, and that u is aC^ local isometry for 
the metrics induced by g on Bud (M U M') and /C fl F (proposition 4.2). Finally, we observe that, for 

r e % n A?", 

a;(^'(r)) = a;(r) = /(a(r)). 
The first equality is due to the fact that 'l'(r) and r are on the same orbit ; the second one is an immediate 
consequence of the definition of the Kasner map /. This shows that the last statement of proposition 6.4 
is equivalent to the following statement about the Kasner map : there exists a constant u > 1 such that, 
for every p G C and every v G TpJC, one has ||D/(p).'y||g > 

This last statement is an immediate consequence of the elementary properties of the Kasner map, and 
of our choice of the riemannian metric g. Indeed, the riemannian metric g was chosen so that it induces 
the same metric on /C n F as the euchdean metric h = dSf + dS^ + + dNf + dN^ + diV| (see (14) 
at the end of section 3). And, since C is a compact subset of the Kasner circle /C which does not contain 
any of the three special points Ti, T2, Ta, there exists a constant uc > 1 such that, for every q £ C and 
every v G TpK., we have \\Df{p).v\\h > i^c\\v\\, where || • \\h denotes the metric induced on /C by the 
euclidean metric h. □ 

Let U := $"^(V). Clearly, is a neighborhood of C n (M U M') in M U M'. Combining proposi- 
tions 6.3 and 6.4, one immediately gets : 

Proposition 6.5. The "Poincare map" Q is well-defined and on U. The compact set C n (M U M') 
is a hyperbolic set for Q. More precisely, there exists a constant v G (0, 1) such that, for every q G 
c n (M U M'), 

• d@{q).v = for every v G F^, 

• de{q).F^ = F^^^y and \\{dQ{q))-\v\\g < v\\v\\g for every v e F^{q). 

Proposition 6.5 shows that the map : M U M' M U M' and the set CnM satisfy the hypotheses 
of the stable manifold theorem 6.2 (for any contraction rate /j, > 0). This shows the existence of local 
stable manifold, with respect to the map 0, for the points of C fl (M U M') : 

Theorem 6.6. For every e small enough, for every q £ C (1 (M U M'), the set 

W^{q) := {r G M U M' I J/5fg(e"(r), G"(g)) < s for every n > 0} 

is a -embedded disc of dimension 2 in M L) M', tangent to F^ at q, depending continuously on q 
(in the topology). For every constant p, > 0, there exists another constant K such that, for every 
g G C n (M U M'), for every r G (g) and every n > 

distg{@'^{r),Q''{q)) < Kp"". 

7. Stable manifolds for the flow 
We are left to prove that theorem 6.6 implies our main theorem 1 .3. 

Consider a point g G C. Then heteroclinic orbit Og j^^g^ intersects the "Poincare section" M U M' 
at one and only one point, that we denote by q. Note that q e C Ci {M U M'). The set VF/ (6, q) 
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defined in the statement of theorem 6.6 is a C^-embedded two-dimensional disc in the three-dimensional 
hypersurface with boundary M U M'. This disc is tangent to F| at q. Since the two-dimensional 
submanifold Byn^ U is not tangent to F| at q, this implies that W^{Q, q) n Bix contains a C^- 
embedded two-dimensional disc in M U M'. Moreover, this disc depends continuously on q. Consider a 
point r e W^{e, q) n i^ix- 



Claim 1. There is an increasing sequence of times {tn)n>0 such that distg{X "(r), f^{q)) — > 0. 

n— >-oo 

According to theorem 6.6, we have 

(31) distg(e"(r),e"(g)) — > 0. 



n— ^-oo 



Together with the continuity of the flow of X, this shows the existence of a increasing sequence of real 
numbers (r„)„>o such that 

(32) dist3(X-"(G"(r)),o;(e'^(g))) 0. 

n— >-oo 

Since r G Bi^, there exists an increasing sequence of times {On)n>o such that, for every n > 0, 

(33) e"(r) = X^"(r). 
Since q G iSjj, we have, for every n > 0, 

(34) uj{Q^m = rHq)) = r^Hq). 

For every n > 0, let tn+i ■= Tn + On- Then (tn)n>o is an increasing sequence, and 

(35) distg(X*"+i(r),r+^(g)) ^ 0. 

n— ^-oo 

Claim 2. The Hausdorjf distance between the piece of orbit {X*(r) ; tn < t < in+i} cind the hetero- 
clinic orbit Ojn^g^ j-n+i^g) tends to when n goes to +oo. 

To prove this claim, we decompose the piece of orbit {X^{r) ; tn <t < tn+i} into three sub-pieces: 

• First, the piece of orbit going from X*"(ro) to $(0"(r)), contained in V. Formula (16), to- 
gether with (35) shows that, for n large, this piece of orbit is close to the heteroclinic orbit 
Ofn{q),fn+^{q)- In particular, for n large, the point $(6"(r)) is close to the heteroclinic orbit 

• Then, a piece of orbit going from $(0"(r)) to G""'"^(r), contained in B \ V. For n 
large, this piece of orbit is close to the heteroclinic orbit C'jn(g) jn+i(g). Indeed, for n 
large, the point ^>(0"(ro)) is close to the heteroclinic orbit Ojn(g) jn+i(g-), and if we write 
e"+i(r) = X*(*(®"W))(^>(e"(r))), then f($(e"(r))) depends continuously on $(e"(r)) 
(proposition 6.4) and thus is uniformly bounded. 

• Finally, a piece of orbit going from (r) to X^"+^ (r), contained in V. Formula (16) together 
with (35) show that, for n large, this piece of orbit is close to the heteroclinic orbit O 

Now, we fix 77 > 0, and we set 

(36) ^fx(^):= U x'{w^{Q,Q)nBix)- 

—r)<t<'q 
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Since W^{@,q)) fi Bjx is a C^-embedded 2-dimensional disc in M, and since the orbits of X are 
transversal to M, we get that Dj^{q) is a C^-embedded 3-dimensional disc in Bjx which depends 

continuously on q. Moreover, the orbit of every point of is the orbit of a point of W^{&, g)nHix- 

So, claim 1 and 2 are still satisfied if we replace r € W^{Q, g) by r G D^lX{q). 

To complete the proof of our main theorem 1.3, it remains to show that the disc D^^^q) depends 
continuously on q. This follows from the continuity of the maps q i-)- W^{Q,q) (theorem 6.6) and 
q^ q (proposition 4.2). 

8. Existence of closed /-invariant aperiodic subsets of the Kasner circle 

The purpose of this section is to prove proposition 1.5. This proposition should be quite obvious for 
people with some culture in dynamical systems. Indeed, the Kasner map / : /C — ^ /C is a degree —2 
map of the circle /C. This implies the existence of a continuous degree 1 map ^ : /C — > M/Z such that 
(f) o f = m o (j) where m : M/Z R/Z is defined by m{6) = —29. Moreover, the Kasner map / 
is expansive : the norm of the derivative of / (calculated with respect to the metric induced on IC by 
h = (iSf + d'El + + diVf + + clN^) is strictly bigger than 1, except at the three special points 
Ti, T2, T3 (where it is equal to 1). This implies that the map : /C — t- M/Z defined above is one-to-one, 
that / : /C — 7- /C is topologically conjugated to the map 6 ^ —29. Finally, it is well-known by experts 
that, for > 2, the union of all compact subsets of M/Z which are aperiodic for the map 9 ^ k9 
is dense in M/Z. We now give a more detailed proof of the proposition for the readers who may not 
necessarily be familiar with low-dimensional dynamics. 

Proof. Recall that Ti,T2,Ts are the three Taub points on the Kasner circle /C. Let Ii, 12,13 be the 
closures of the three connected components of /C \ {Ti, T3}, the notations being chosen so that Ti is 
not one end of /i, T2 is not one end of I2, and T3 is not an end of /a. We denote by dist the distance on 
/C induced by ft. = dSf + dS^ + dT^j + dN^ + dN^ + dN^ ; the norm of the derivative of / will be 
calculated with respect to this distance. 

We consider the set S := {1,2, 3}^ endowed with the product topology, and the shift map a : S — )• S 
defined by a{ao, ai, 02, ■■ ■)i = (ai, 02, 03, . . . ) (in other words, if a = (ai)igp} G S, then {a{a))i = 
Oj+i). Let So be the subset of S defined as follows : 

So = {a = (ai)igp} G S such that a^+i 7^ ai for every i}. 

Note that So is cr-invariant. We will construct a continuous "almost one-to-one" map : Sq ^ /C such 
that (f)o a = f o (p. 

Claim. For each sequence a := (ai)igN in Sq, there exists a unique point x & IC such that /'(x) G /„. 
for every i >Q. 

In order to prove the existence of x, one just needs to notice that the image under / of each of the 
intervals Ii, I2, 13 is the union of the two other intervals. This implies that the intersection Hi^o f~^{Iai) 
is non-empty for every N, and therefore, that the intersection HieN f^^i^a^) is non-empty. The existence 
of X follows. In order to prove the uniqueness of x, observe that: for every £ > 0, there exists a(£) > 
such that II \Df{x) || |/i > 1 + a{e) for every x e K. such that dist(a;, Tj) > e for i = 1, 2, 3. Hence, if 
X ^ x' were two points such that f^{x) G and G for every i > 0, then one would have 
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dist(/*(a;), p{x')) oo ; this is absurd since the lengths of /i, I2 and I3 are finite. Hence there is at 
most one point x in Dign f~^{Iai)- This completes the proof of the claim. 

Now, we consider the map : Sq — )■ /C which maps a sequence a := (aj)jgN to the unique point 
X e }C such that /*(a;) G Iq. for every i > 0. This map obviously satisfies (poa = focp. It is continuous 
(this is an immediate consequence of the continuity of /) and onto (because the image under / of one of 
the intervals /i, /2, 13 is contained in the union of the two others). It is not one-to-one. For example, the 
Taub point T3 has two pre-images under (p : the sequences (1, 2, 1,2,1,2,...) and (2, 1, 2, 1,2,1,... ). 
More generally, every point x e IC such that f''° (x) is a Taub point for some integer io > (and thus 
P (x) = /*° (x) for every i > zq) has two pre-images under (f>, and these two pre-images are preperiodic 
for a. This is the only lack of injectivity of : if x € /C has not a single pre-image under 0, then there 
exists io > such that /*o (x) is a Taub point (this follows from the fact that the intersection between 
two of the three intervals hjh, h is reduced to a Taub point). 

It follows that the image under of a closed cr-invariant aperiodic subset of is a closed /-invariant 
aperiodic subset of /C. So we are left to prove that the union of all the closed cr-invariant and aperiodic 
subsets of TiA is in dense in S^. 

A element a = (aj)j>o of Sq is said to be square-free if it does not contain the same word re- 
peated twice : for every io > and every £ > 0, the word . . . aig+^-i is different from the word 
ajp+^ . . . ai„+2^ 1- It is well-known that there exist square-free elements in Sq (such an element may be 
easily deduced from the well-known Prouhet-Thue -Morse sequence, see for example [1, corollary 1]). 
Now let a = {ai)i>o € Sq be square-free, then the a-orbit of a does not accumulates on any periodic 
(T-orbit, hence the closure of the a-orbit of a is a closed, positively invariant aperiodic subset of Sq. 
Moreover the same is true, if one replaces a by a sequence a' G Sq which has the same tail as ct (i.e. 
there exists io such that a' = a-i for i > io). The set of all sequences a' which has the same tail has a is 
obviously dense in Sq. Hence the union of all closed cr-invariant aperiodic subsets of Eq is dense in Sq- 
As explained above, the proposition follows. □ 
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